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Abstract 



Doppler shift measurements of the Na Di absorption line have revealed solar oscilla- 
tions in a new regime of frequency and wavenumber. Oscillations of vertical velocities 
in the temperature minimum and low chromosphere of the Sun are observed with fre- 
quencies ranging up to 9.5 mHz. There is no evidence for chromospheric modes of 
3 minute period. This indicates that the chromosphere does not form a good cavity 
for acoustic waves. 

The fundamental-modes appear with wa\^enumbers up to 5.57 Mm“^ (equivalent 
spherical harmonic degree, 3S77). The frequencies lie below the predicted values at 
wavenumbers above 1 Mm”b The values are in agreement with previous measure- 
ments that exist for wavenumbers up to 2.67 Mm“^ 

Spatial maps of velocity power show that high wavenumber oscillations are sup- 
pressed in active regions. The shape of the power depression indicates that wave 
motion is affected in the layer of atmosphere where the measurement is made. The 
/-modes are suppressed in the same way as p-modes, indicating that the mechanism 
for wave suppression affects velocity fluctuations. Mode frequencies are not affected 
by the magnetic fields by more than 50 /iHz, the precision of the measurement. 



IV 



Acknowledgments 



In the course of my graduate student career at Stanford, I have made many friends, 
and benefited from the knowledge and experience of many people. My advisor, Phil 
Scherrer, has always made himself available for discussions, and shared his insight in 
the field. He has also provided guidance for this dissertation from beginning to end. 
I am grateful to him for encouraging nty participation in workshops and conferences 
while, at the same time, giving me the freedom to explore my interests outside the field 
of solar physics. Dave Burke, my former advisor at the Stanford Linear Accelerator 
Center, deserves thanks as well, for giving me a solid grounding in experimental 
physics. 

Ted Tarbell and Alan Title at the Lockheed Palo Alto Research Laboratories 
(LPARL) have both taken a personal interest in this project. They provided the 
instrument and expertise to make the observations, as well as many ideas and topics 
for research. They also offered many useful comments on this manuscript. All of 
the data analysis was done at LPARL, with support from Kermit Smith and Zoe 
Frank. I am especially grateful to Dick Shine for his help and consultation in software 
development. 

I would like to thank Douglas Gough, Pawan Kumar, and Ed Lu for illuminating 
discussions that have greatly increased my understanding of helioseismology theories. 
Thanks also to Ellen Zweibel, Neil Hurlburt, and Colin Rosenthal for interesting 
discussions of the results, to Tim Brown, Goran Scharmer, and Louis Strous for 
their suggestions concerning data analysis, to Paco .Armas and Rolf Kevar for their 
support at the Swedish Solar Observatoiy, and to Bob Wagoner for his suggestions 
for improving this manuscript. Margie Stehle, Gail Ragsdale, and Marcia Keating 



have good-luimorecll}- and patiently guided me through the mountains of paperwork 
and miles of redtape of Stanford administration. 

I have enjoyed my yc'ars at Stanford, and learned a great deal thanks to the 
friendship and support of some wonderful people. My classmates and peers in physics, 
Anna Green, Bruce King, Robin Sz Pete Wright, Lee Anne Kowalski, and Steve 
Pollock helped me keep my sanity through a year of problem set hell, and The Qual 
Interaction with Columbaens, Synergites, and assorted new-age hippies has re- 
shaped my views and changed my long-term goals. Suzanne Jones, Sally Otto, Matt 
Price, Neil Switz, and my sister, Karen, are all awesome, and I hope to be just like 
them when I grow up (without the long hair, though). Darlene Zanker, Greg Hollon, 
and Lise Schickel were part of a most memorable summer class at the Hopkins Marine 
Station, as was our instructor, Jim Wataiiabe. Leslie Roldan, Suttirat-Anne Larlarb, 
Sunit Mahajan, and Ann Fridlind from (.’olumbae, and Ethan Pride, Cecil Scheib, 
and Myia Johnson from Synerg\* ha\'e all been willing to listen to my grumbling about 
'Writing and, to provide me with distractions and diversions. 

Finall}^ I wish to tliank Peter \htousek, Carla D’Antonio, and Jim Raich for 
advice and help getting me on my feet in a new field. 

This work was funded b\* NAS.A contracts for SOUP (NAS S-32S05), OSL-CIP 
(NAS 5-26S13), and SOl/MDl (NAS 5-303S6), by NASA grant NGR-559, NSF grant 
ATM 90-22249, ONR grant N00014-S9-J-1024, and by Lockheed Independent Re- 
search Funds. 



\'i 



Contents 



Abstract iv 

Acknowledgments v 

1 Introduction 1 

1.1 Wave Properties 2 

1.1.1 Acoustic Waves 2 

1.1.2 Surface Waves 4 

1.2 High Frequency .Modes 6 

1.3 High Wavenumber .Modes 8 

1.4 Height of Observations 9 

2 Measurement Techniques 11 

2.1 Three-Point Measurements of Line Parameters 12 

2.1.1 Algebraic Solution 13 

2.1.2 Error Estimates 13 

2.2 Four-Point Measurements of Line Parameters 14 

2.2.1 Implementing the .Algorithm 17 

2.2.2 Spectral Simulation 20 

2.2.3 Sources of Systematic Error 21 

2.2.4 Statistical Error 25 

2.2.5 Discussion of the Method 29 

2.3 Magnetograms 30 

2.3.1 Inverse Zeeman Effect 30 

vii 



2.3.2 Loiigit iiclinal Field Measurement 30 

2.4 Transform Methods 31 

2.4.1 Practical Considerations 33 

2.4.2 .Apodization 34 

2.4.3 Normalization 35 

2.4.4 Statistical Error Estimates 35 

3 Telescope and Instrumentation 36 

3.1 Telescope 36 

3.2 Postfocus Instrumentation 36 

3.2.1 Tunable Filter 36 

3.2.2 CCD Camera 42 

3.2.3 Image Stabilization System 43 

3.3 Optical Layout 44 

3.4 Modulation Transfer Function of the System 48 

3.4.1 Entrance Pupil 49 

3.4.2 CCD Camera 49 

4 Observations 51 

4.1 Observing Procedure 51 

4.2 Calibration .Method 52 

4.3 Solar Observations 55 

5 Data Reduction 61 

5.1 Flat Fielding 62 

5.2 CCD Flaw Corriction 63 

5.3 Image Rotation 63 

5.4 Atmospheric Distortion 63 

5.4.1 MTF Normalization 64 

5.4.2 Destretching 66 

5.5 Dopplergrams and Intensit ies 67 

5.6 Image Drift 68 

viii 



5.7 Angular Scale Delenninatioa 



68 



6 Measurements from Power Spectra 71 

6.1 Collapsing 3-cl Power Spectra 74 

6.1.1 Distance Scale Determination 74 

6.1.2 Azimuthcil Projection 76 

6.2 Parametrization and Least-Squares Fitting 78 

6.2.1 Fitting Procedure 78 

6.2.2 Levenberg-.Marquardt .Method for Nonlinear Regression .... 80 

6.3 Results of Fits 82 

6.3.1 Background 82 

6.3.2 Mode Frequencies 83 

6.3.3 Mode Widths 87 

6.3.4 Mode Power 88 

6.3.5 Comparison of Magnetic and Nonmagnetic Regions 91 

7 Spatial Distribution of Power 94 

7.1 Mode Masks 94 

7.2 Results 95 



8 Conclusions and Discussion 102 

8.1 Surfaces and Ca\ities 102 

8.2 Waves in Magnetic Fields 103 

8.3 Directions for Future Work 104 

A Frequency Tables 106 

Bibliography 112 



IX 



List of Tables 

4.1 Filter Parameters for SEISX.AFE Sequence 56 

4.2 La Palma Observing Runs using SEISN.AFE Secjuence 57 

4.3 Filter Parameters for .XAF.XS'f Sequence 57 

A.l Mode Frequencic's from Dopplergrams 106 

A. 2 /-mode Frec(uencies from Intensities 110 






List of Figures 



1.1 Acoustic ciitoir^v vs. /,• 3 

1.2 Propagatiou of an acoustic wave in a plane-parallel atmosphere. From 

Gough (1983) 4 

1.3 Temperature and density of the solar atmosphere. From Avrett (1992). 5 

1.4 Mode frequencies vs. wavenumber with a. chromospheric cavitj' present. 

From Ulrich lU Hhodes (1977) 7 

1.5 Profiles and contribution functions for the Na D lines. From Schleicher 



(1976) 10 

2.1 Measurements of the .\'a Di Hue 12 

2.2 Spectral profiles 15 

2.3 Ratio o vs. velo( ity 16 

2.4 Ratio j3 vs. velocity 17 

2.5 Transmission product for Lyot and Michelson filters vs. wavelength . 22 

2.6 Transmission ratio //adA/ / f\d\ as a function of fwhin of the Lyot 

filter 22 

2.7 Systematic error t,, vs. velocity for line-depth variations 23 

2.S Systematic error t^, vs. velocity for Lyot drifts 24 

2.9 Systematic error c^. vs. velocity for line-width increments: 0, 10, 20, 

30% (o-onlj' met hod) 24 

2.10 Systematic error t„ vs. velocity for line-width increments: 0, 10, 20, 

30% (o, d method) 25 

2.11 Geometry of modeled line asymmetiw 26 



.\'i 



2.12 Depth of measured velocity vs. x'elocity for bisector shape hwhm: 

— 7.5 ,7.5 mA (a-oiily method) 26 

2.13 Depth of measured \-elocity vs. velocity for bisector shape hwhm: 

— 7.5 7.5 luA (a, 3 uu'thod) 27 

2.14 Statistical error (7^, vs. velocity for different Hue-widths (a-only method). 28 

2.15 Statistical error vs. x^elocit}’ for different line-widths (a, ^ method). 28 

2.16 Measurement of Zeeman splitting 31 

3.1 Schematic diagram of the Swedish 50 cm Vacuum Solar Telescope. 

From Scharmer tt al. (1085) 37 

3.2 Swedish Solar Ol>ser\'atory. La Palma (photograph) 38 

3.3 Pulse propagation through a simple element. From Title & Rosenberg 

(1981) 39 

3.4 Lyot filter with four elements. From Title Sz Rosenberg (1981) 39 

3.5 Transmission vs. frecjiieucy for a L}ot filter. From Title Sz Rosenberg 

(1981) 40 

3.6 Birefriugent tuning ehnnenl. From Title Sz Rosenberg (1981) 41 

3.7 Transmission profile' of the SOFP tunable filter 43 

3.8 Relatix'e transmission of the SOFP tunable filter 44 

3.9 Optical bench se'lup (photograph) 46 

3.10 Optical layout at Swc'dish Solar Observatory 47 

3.11 Timing offset for images 48 

3.12 Sc[uared MTF of teh'seope and C'CD detector 50 

4.1 Line scan of Xa \)\ Line using the CCD camera 53 

4.2 Full-disk magnet ogra in from KPXO 58 

4.3 Continuum iuteusily image' of ARGISI 59 

4.4 La Palma magiu'te>gram of ARG18I 59 

4.5 H a image of ARG181 60 

5.1 Spatial frequeucx' log-spectra of Na Di —40 mA filtergrams 65 

5.2 Image drift during obserx’ing run 69 



6.1 Magiietogram of ARG181 showing siibfielcls for study 72 

6.2 Surfaces of power in (A\/^) space 73 

6.3 Constant frequency slice of the 3-d power spectrum 75 

6.4 Scaled velocity-power spectrum. fMom Fernandes et al. (1992) .... 77 

6.5 Constant wa\*ennmber slices of the 3-d power spectrum 79 

6.6 Power law e.xpoiient for background fit 82 

6.7 Background \elo(dty power integrated o\er frequency. From Fernandes 

e/a/. (1992) 83 

6.8 Calculated Doppler \*elocity vs. lieight. From Nesis &: Mat tig (1989) . 84 

6.9 Measured mode frec[iK'ncies (from velocity spectra): i/ vs. log^o * 85 

6.10 Comparison of measured p-mode frequencies with Libbrecht et al. (1990) 86 

6.11 Comparison of /-mode frequencies measured from velocity and intensity 87 

6.12 Comparison of measured /-mode frequencies with Libbrecht et ai (1990) 88 

6.13 Comparison of measured /-mode frequencies with theory 89 

6.14 Frequency and ridge width vs. k 90 

6.15 Velocity power \ s. mode frequency. From Fernandes et al. (1992). . . 91 

6.16 Mode frecjiiency differences for canopy and quiet subfields 92 

6.17 Mode frequency differences for sunspot and quiet subfields 93 

7.1 Masks for power distribution maps 96 

7.2 Comparison of power maps with magnetic fields 97 

7.3 Power distribution for selected masks 98 

7.4 RMS velocity vs. k 99 

7.5 Slice of /-mode map 100 

7.6 Comparison of /-mode slice with magnetic field 101 



xiii 



I 



! 



i 



C 



Chapter 1 
Introduction 



Solar oscillations were first observed by Leighton (1960), and were explained by Ulrich 
(1970) and Leibacher & Slein ( 1971) as being the superposition of many global stand- 
ing waves. They ha\'e subsecjnently become an important tool for the determination 
solar structure. Measurements from helioseismology include the depth of the con- 
vection zone (Gough 1976), helium abundance in the convection zone (Christensen- 
Dalsgaard & Gough 19S0), and internal rotation rate (Duvall et al. 1984; and others). 

In the past, obser\’ational limitations have restricted helioseismology to long wave- 
length oscillations (~ 4 Mm). New instruments and technicjues are beginning to make 
shorter wavelengths accessible to study. This opens the door to the exploration of 
structure at much smaller spatial scales. Stratification of the atmosphere, and vertical 
structure of magnetic regions are becoming active areas of investigation. 

The term solar osciUations conv’entionall}^ refers to two types of waves. Because 
they exist in a resonant cavit}', these waves are also called modes (eigenmodes of 
the cavity). Compressional or acoustic waves (p-modes) are longitudinal waves with 
pressure as the restoring force. The}^ have at least one nodal surface in the radial 
direction. Gravity waves ha\e grax iU' as the restoring force. A special case of these 
are the surface waves. Bc’cause they ha\'e no radial nodes, surface waves are called 
fundamental modes (/-modes). These transverse waves propagate at the interface 
between two media. The remainder of tlie gravity waves are called p-modes. To date, 
p-modes have not been unambiguously detected, and they will not be considered here. 
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1.1 Wave Properties 

1.1.1 Acoustic Waves 

Sound waves can propagate vertically at frequencies greater than the acoustic cutoff 
or critical cutoff frequency At lower frequency, a disturbance moves the whole 
atmosphere up and down in jihase. Hence, waves do not propagate. At low horizontal 
wavenumber k (^k Hie acoustic cutoff can be expressed in terms of the adia- 

batic sound speed c and ihe density scale height H. For an isothermal atmosphere, 
the scale height is given by 

ll = - (1.1) 

19 

and, the acoustic cutoff, liy 

2H 

where g is the acceleration due to gravity; and, 7 is the adiabatic index (Lamb 
1932). Taking c = 7.1 kin/s. g — 27S m/s^, and 7 = | in the photosphere, we find 
H = 110 km and = 5.2 mHz for h 4.5 Mm“H 

In general, the acoustic cutoff depends on wavenumber. The relationship is given 
by (Kumar iC Goldreich 19S9) 

= ( 277 )' ! Ill + + y^(l + (2A-//)T-§(2t//)4 (1.3) 

Figure 1.1 shows a function of k. 

Oscillations with wavelength greater than a few Mm propagate around the sun 
without excessive dissipation. A periodic boundary condition is, therefore, imposed 
on these waves. The wa\'e equation is separable into radial and angular parts and has 
eigensolutions of the form, 

vI/„;„d'-.Co,/) = R,Ar)Yr{0,<l^)e-‘^^ (1.4) 

n = 1,2,3.... / = 0,1,2,... m = —/,—/+ 1 , ..., / 



Short wavelength oscillations are dissipated before traveling around the sun. More- 
over, at their small spatial scales, the inecliiim can be approximated as plane-parallel. 
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0 2 4 6 8 10 

k (Mm“^) 

Figure 1.1: Acoustic cutoff ujc vs. k 



Plane waves, therefore, pro\*icle a more natural basis set for the horizontal eigenfunc- 
tions. 

= (1.5) 

The horizontal wavenumber h = |A-| can be related to the spherical harmonic degree 
/ as follows. 

k = ^ ^ (1.6) 

^sun 

Ray paths for a p-mode are diagrammed in figure 1.2. As waves travel downward, 
they encounter higher temperatures and an increasing sound speed. This causes them 
to refract back upward wliere they encounter increasing values of Uc- When the wave 
can no longer propagate (^’ < u.v-). it is reflected back down. The resonance condition 
is met when crests of downgoing waves are superposed on crests of refracted waves 
returning to the surface. For waves trapped in an adiabatic and polytropic layer 
(polytropic inde.x = 4), this condition gives a dispersion relation of approximate 
form, 
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Figure 1.2: Propagation of an acoustic wave in a plane-parallel atmosphere. The wave 
is refracted as it propagates downward, and reflected at the surface. The maximum 
depth reached by the wa\*e is called the lower turning point. From Gough (1983). 



w-'=(| + l)</A; u = 1,2,3,... 
where g is the acceleration due to gravity at the surface. 



(1.7) 



This dispersion relation can be measured directly from the ridges of power in a 
spectral decomposition (a plot of power vs. and k) of solar velocities or intensities. 
The exact relation depends on the variation of sound speed with depth. 

The lower turning point d of the waves is approximately given by. 



d = 



3 ( 7 ! + 2 ) 

ik 



( 1 . 8 ) 



The excitation source of the p-modes is thought to be turbulent convection at 
the top of the convection zone. The efficiency of acoustic emission from convection 
is proportional to approximately the eighth power of the Mach number (Goldreich 
& Kumar 19SS). Hence, only a few of the granules on the high velocity tail of the 
distribution are expected to contribute to the acoustic noise driving the p-modes 
(Brown 1991). 



1.1.2 Surface Waves 

Surface waves can propagate on the interface between convectively stable media (e.g., 
air over water). On the Sun. the transition region between chromosphere and corona 
provides such a boundary (.Jones 1969). The model of Avrett (1992) describes the 
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Figure 1.3: Temperature and density of the solar atmosphere. From Avrett (1992). 

temperature and density discontinuities shown in figure 1.3. For purely vertical den- 
sity stratification of a compressible fluid, surface waves obey the dispersion relation 
(Whitney 1958), 

= gk (1.9) 

where g is the local acceleration due to gravity. 

The simple form of this relation has prompted the remark that, “There is no 
interesting physics in the /-modes.'’ It wasn’t until /-modes were found not to obey the 
above relation at high wavenumber that this conclusion was reexamined. Libbrecht, 
Woodard & Kaufman (1990) found /-mode frequencies to lie below the predicted 
values for wavenumbers above 1.4 .\Im~^. .At least three possible explanations exist, 
though numerical predictions have not yet been offered. 

Campbell & Roberts (1989) and Evans & Roberts (1990, 1991) predict that a 
horizontal magnetic field in the chromosphere would increase the frequencies of the 
/-modes by increasing the restoring force on the surface waves. The sign of the 
frequency shift may reverse at high wavenumber, however, in a more realistic model 
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of the chromosphere with an appropriate magnetic field structure (Roberts 1991). 

Rosenthal (1992) finds that evanescent wav'es dissipated in the corona would cause 
power to leak out of the /-modes. This would lower the frequency of the modes at a 
given wavenumber by an amount dependent on the temperature difference between 
chromosphere and corona. 

Finally, structures in the atmosphere, such as flux tubes and convection cells, 
violate the assumption of horizontal homogeneity. It is not known whether this should 
increase or decrease the mode frequencies (Gough 1992). 



1.2 High Frequency Modes 

For the purpose of cliscus.sion here, the high frequency regime includes modes with 
io > where = 5.2 mHz. This value of uJc is the photospheric value for low 

wavenumbers. Because wav^es with can propagate, they are not trapped by 

the photosphere. So, in the model illustrated by figure 1.2, there should be no upper 
reflecting surface at frequencies > 5.2 mHz. Removal of the top surface of the cavity 
eliminates the possibility of resonance, and hence, the mode structure. 

It was a surprise, thendbre, when the /;-mode ridges were found to continue above 
the acoustic cutoff (Jeffries et al. 19SS; Libbrecht 19SS; Duvall et aL 1991). Two expla- 
nations have been put forward for the origin of the high frequency ridges. Balmforth 
h Gough (1990) propose that strong reflection from the base of the corona defines a 
chromospheric cavity which couples to the subphotospheric cavity. Waves with fre- 
quency high enough to pass through the photosphere would still be trapped by the 
corona. 

The presence of a chromospheric cavity has further consequences, however. The 
chromosphere is of order one wavelength in thickness for 5-minute modes; and waves 
of appropriate wavelengtli A,, will resonate in this cavity. The p-modes propagate 
almost vertically in the atmosphere (angle of incidence 0). Thus, A,. = cjvr^ 
independent of horizontal wavenumber. The chromospheric modes should, therefore, 
have a dispersion relation. 



; = const. 



( 1 . 10 ) 
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Figure 1.4: Mode frequencies vs. wavenumber with a chromospheric cavity present. 
From Ulrich Rhodes (1977). 

Ulrich Sz Rhodes (1977) found that the p-mode ridges would merge with a set of 
chromospheric modes at about 5.5 mHz. See figure 1.4. The p-modes would continue 
on abov'e the 5.5 mHz band with an '‘avoided crossing” between the ridges above and 
below the chromospheric resonance. 

Dame, Gouttebroze h Malherbe (19S4) claim to have seen the chromospheric 
modes in Ca II K-line observations. Their power spectra are very noisy, however, and 
it is difficult to discern any definite ridges of power. Deming et al. (1986) find a peak 
in the one-dimensional power spectrum (power versus uj) at 4.3 mHz. Its significance 
is difficult to assess, though, without knowledge of the wavenumber dependence of the 
power. Woodard Iz Libbrecht (1991) do not find this peak in similar one-dimensional 
power spectra. 

The other explanation for the high frequency ridges was offered by Kumar et al. 
(1990). They propose that waves are emitted isotropically from a source near the top 
of the convection zone. W aves traveling directly toward an overlying observed surface 
interfere (in Fourier space) with wa\’es refracted upward from below to produce “mock 
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modes.” It must be strcssc'd that this interference is not physical. In fact, the directly 
emerging \va\^e trains and refracted wave trains may arrive at the surface at different 
times and in different place.s. The basis functions, used to represent the waves, can 
only distinguish time, modulo 2i: juj, and position, modulo jk. Frequencies of the 
mock modes can be used to determine the depth of the acoustic source layer (Kumar 
k Lu 1991). 

The two models for high frequency waves can be distinguished by observations 
that resolve the ridge structure in the 5.5 mHz frequency range. 



1.3 High Wavenumber Modes 

Most of the literature in helioseismology refers to the (spherical harmonic) degree 
/ of a mode rather than the wavenumber k. By convention, low degree refers to 
0 < / < 4, intermediate degree, to 4 < / < 100, and high degree to 100 < /. The 
oscillations presented in this dissertation have 500 < / < 3877. By analogy with 
the electromagnetic spectnim. the\' should, perhaps, be called very high or ultrahigh 
degree, but I will stick to the term, high wavenumber. 

Previous obserx'ations of the /-modes exist up to wavenumbers of 2.67 Mm“^ 
(Libbrecht, Woodard k Kaufman 1990). As mentioned above, the frequencies do 
not agree with equation 1.9 at wa\’enumbers above 1.4 Mm"^. The variation of g 
with height accounts for onl\' a 0.5% modification to equation 1.9. This is an order 
of magnitude less than the obserx’ed effect. It would be desirable to confirm the 
observational result and determine the trend at higher wavenumbers. 

The size of structures that can be resolved by oscillations is roughly equal to 
the wavelength of the modes. For actix’e regions, spatial scales range from 200 km 
(flux tubes) to 25 Mm (large sunspots). The corresponding wavenumbers range from 
0.25 to 30 Mm*”F Althougli the high end of this range is not attainable with any 
current instrument, there is still much room for investigation. 

Woods k Cram (1981) have reported a suppression of oscillation amplitudes in 
plage as compared to qui(4 regions. Tarbell et aL (1988) and Title et aL (1992) find a 
suppression of both conwcti\*e motions and wave amplitudes. The suppression of the 
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oscillations coincides with the boundary of the photospheric magnetic field to within 
one wavelength. 

Braun, Duvall Sz LaBonte (19S7, 19SS) and Braun et al. (1992b) have reported 
absorption of /;-mode power by sunspots. The mechanism for this is not understood. 
Resonant absorption (Hollweg 19SS; Chitre k. Davila 1991), scattering to short wave- 
length (Bogdan k Zweibcl 19S7), and enhanced viscous dissipation (LaBonte k Ryu- 
tova 1992) are among the explanations offered. It is by no means clear whether 
only propagating acoustic waves below the photosphere are absorbed, or whether 
evanescent waves are also affected. A spatial map of oscillation amplitudes above 
the photosphere, where the magnetic fields are different, would help untangle the two 
effects. 

Global oscillations ha\*e been found to have frequency shifts correlated with solar 
activity (Palle, Regulo k Roca Cortes 1989; Woodard et al. 1991). This may be due 
to a change in the shape' of the ca\dty surface with sunspots present. Waves of a 
given frequency may have different wavenumbers within magnetic and nonmagnetic 
regions. Global measurements would then find a weighted average of the two disper- 
sion relations. Comparison of high resolution observations in quiet sun and active 
regions could resolve this issue by direct measurement of dispersion relations in the 
two types of regions. 



1.4 Height of Observations 

Oscillations can be obsei ved as fluctuations in intensity and line-of-sight velocity. 
Both p-modes and /-modes produce motions of gas particles, and cause Doppler 
shifts of absorption lines. Pressure fluctuations due to p-modes produce tempera- 
ture/intensity fluctuations. The /-modes also produce intensity fluctuations, though 
the reason is not well understood. It may be that the ~ 100 km vertical displacements 
of parcels of gas affect th(' process of absorption line formation. 

Different parts of an absorption line are formed in different layers of the atmo- 
sphere. Hence, the height at which the velocities and intensities are determined 
depends on the positions in the absorption line at which measurements are made. In 
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Figure 1.5: Na D line profiles on left hand scale. Contribution functions for emergent 
intensity on right hand scale. Thick lines show the maximum of the contribution 
function at each wavelength. Thin lines show 1/3 of the maximum of the contribution 
function. From Schleicher (1976). 

the core of deep lines such as Na Di, very little continuum light emerges from the 
photosphere. The height of a measurement, therefore, depends in the contribution 
function for the emergent line emission (Gurtovenko, Sheminova & Sarychev 1991; 
Schleicher 1992). This function is shown in figure 1.5. 

The original intention of the observations described in chapter 4 was to make one 
set of measurements in the core and one in the wings of the Na Di line. Problems 
described in section 4.1 resulted in a measurement weighted mostly in the core (600- 
850 km), but with some contribution from below the temperature minimum (150- 
500 km). 



Chapter 2 

Measurement Techniques 



Producing power spectra from uncalibrated images involves several stages of data 
analysis and reduction. 1 his chapter re\dews some general procedures, standard and 
nonstandard, for measuring line-of-sight components of velocity and magnetic field, 
and for determining wave amplitudes and phases. The details of calibration to remove 
artifacts of the instrument and atmosphere are left to chapter 5. 

The line-of-sight velo( ity is found from the Doppler shift of an absorption line. 
The Doppler shift, in turn, is determined b}^ measuring intensity as a function of 
wavelength at a series of fixed-wavelength positions in the absorption line. The more 
points in the measurement, the better the velocity estimate will be. Unfortunately, 
more points require proportionately more time to make the measurements, and geo- 
metrically more time to analyze them. 

In this chapter are descril)ed two new technicpies for making velocity images (dopp- 
lergrams) and line-center intensity measurements using three-point and four-point 
methods. The four-point iiK^tliod was developed for the Solar Oscillations Investi- 
gation / Michelson Doppler (SOl/MDI) experiment for the Solar and Helio- 

spheric Observatory (SOHO) spacecraft. It was not applied to the data used for this 
dissertation. 

Images of line-of-sight magnetic field strength (magnetograms) are an important 
tool in this study. The tried and tested “two-point magnetogram” method is described 
below. 
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Figure 2.1: Measurements of the Na Dj line 

Finally, the Fourier transform is essential to the study of waves propagating in 
a linear medium. Procedures used to extract information about the oscillations are 
described in the last section of this chapter. 

2.1 Three-Point Measurements of Line Parame- 
ters 

The three-point method for finding Doppler shifts and line center intensities assumes 
a parabolic absorption line profile. This assumption is valid when the three intensity 
measurements are positioned relatively close to the line core. Figure 2.1 shows the 
measurements schematically for the Xa Dj (5896 k) line. The Doppler shift and line 
center intensity are the x and y coordinates of the vertex, respectively. 
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2.1.1 Algebraic Solution 



Assuming a parabolic line profile 



t ( A ) — /o + 



A-Aq^ 
const. , 



( 2 . 1 ) 



we can solve for the vert(\\ (Aq, /o) in terms of three intensities Ii, h measured at 
wavelengths Ai^A^.Aa, respect i\’ely. 



^ + ^^ 3^3 

Au — ^ 

C\ + C3 

, ^ , (^iCi - 6303 )^ 

4S,Ss(Ss-SC{c^-Cr) 



where Aj = i(A] + A2), 


^^3 


= |(A 3 + A2), 


Ai = Ai — A2. 


^3 


= A3 — A2, 


CN 

1 

II 


, and C3 


= <^ 3 (/l -I2). 



( 2 . 2 ) 

(2.3) 



2.1.2 Error Estimates 

The errors on the measured ciuantities /i,/ 2,/3 are due to limitations of the instru- 
ment. The transmission of the tunaiile filter \*aries with wavelength as shown in 
figure 3.8. The phase of ('ach Fourier component of this modulation is independent, 
and drifts with time. In addition, the wavelengths tuned by the filter are quantized 
in units of the smallest step size (3.TI mA). These sources of error were used in a 
simulation to estimate errors for the \^elocity and intensity. 

The simulation assumed a parabolic line profile with parameters measured from 
the data. A range of v^elocities of ±2000 m/s with a flat distribution, and a range of 
line-center intensities of ±S9( with a Gaussian distribution were used in the model. 
Intensity '‘measurements** /], T, I 3 were modeled assuming a random jitter of up to 
±1.57 mA (flat distribution) in the location of each measurement, and assuming 
phases for the transmission modulation which varied randomly between sets of 3 
filtergrams, but not within a set. The velocities and intensities were then calculated 
using equations 2.2 and 2.3. Tlu^ errors had approximately Gaussian distributions 
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with standard deviations. 



= 300 m/s, and <7/ = 7.7%. 



(2.4) 



2.2 Four-Point Measurements of Line Parame- 
ters 



An algorithm was cle\’elo])ecl for the SOI/MDI experiment to determine Doppler shifts 
using four measured intensities. The measurements are used to solve a system of si- 
multaneous equations for the four absorption line parameters. Since computing time 
must be kept to an absolute minimum, the nonlinear equations are solved (numeri- 
cally) using a look-up table approach. 

The model included only the Lyot and Michelson filter elements, and the following 
parameters were used: 

• Lyot filter: profile - - Gaussian 

fwhm - 500 mA (see section 2.2.3) 



• Michelson elements (2): profile — cos^{kX) cos’(2A'A) 

fwhm - 105 mA 



/here k = 



fwhm 



cos 



-1 •>-; 



+ \/W)' + - \/lf)' 



Ni I (6768 A) Line: profile 
fwhm 
depth 



Gaussian 
130 mA 
5S^T 



• Measurement Spacing — 85 niA (s\mimetric about nominal line center) 

See figure 2.2 for illustration. 

Four measurements are made at increasingly longer wavelengths across the ab- 
sorption line, to yield intensities Ii through G. The ratios, a and are formed as 
shown: 



CHAPTER 2. MEASUREMENT TECHNIQUES 



15 




-800 -400 0 400 800 
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Figure 2.2: Spectral profiles. 



J = ^ (2.6) 

The quantities a and 3 are independent of CCD dark current and gain variations. 
They are also independent of continuum intensity and line-depth if a suitable choice 
is made for the Lyot filter profile (see section 2.2.3). 

We can solve exactly for the Doppler shift Aq, line depth d, and continuum inten- 
sity Iq if we assume a cosine shaped absorption line with meaisurement spacing equal 
to the hwhm: 



/(A) = 7o 




( 7t(A - Ao) 
\ jwnm 




(2.7) 



Solving, we find, 



Aq = fwhm ■ 



( - — cot * a 
V4 



/o=:^(/i+/3-^^) 



( 2 . 8 ) 

(2.9) 



d = 



/o 



where s = + [I, - 1^)2 



( 2 . 10 ) 
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velocity (m/s) 



Figure 2.3: Ratio a vs. velocity for width increments: 0, 10, 20, 30%. 



Note that a has a one-to-one correspondence with the ratio r = 
which appears in the SOI/.MDI proposal (i.e., r = (a — l)/(a-|-l)). Hence determining 
the velocity from a alone is eciuivalent to determining it from tan“' 

For a Gaussian line-profile. 



r. 



/(A) = lo 



1 — de' 






2<t 



¥ 



( 2 . 11 ) 



where cr = fwhm/2.36 



The ratios q and ^ are plotted as functions of velocity for different line-widths in 
figures 2.3 and 2.4, respectively. It is apparent that the velocity determined from 
a alone would be in error if the calibration table were not made for the ‘‘right” 
line-width. 

This caveat can be removed by determining the velocity v from both a and /?: 
V = u(a,/3). An algorithm for doing this is outlined below. 
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velocity (m/s) 



Figure 2.4: Ratio 3 vs. velocity for width increments: 0, 10, 20, 30%. 

2.2.1 Implementing the Algorithm 

Look-Up Table Approach 

It is possible in principle to compute velocity by using a two-dimensional look-up table 
in a and (3 (i.e., v — v(o,/3)). This, however, would require an inordinate amount 
of memory ((12 bits in q) x (9 bits in 3) = words) to achieve an accuracy of a 
few m/s. A better strategy is to take advantage of the fact that a, alone, provides a 
pretty good estimate of velocity; and write: 

i*(Q,Ai) = i’i(a) -I- V2 (q,/3) (2.12) 

where V 2 is a correction to the velocity t'j and, hence, does not vary with a as quickly 
as does vi. 

This approach can achieve an accuracy similar to the previous approach but with 
much less memory overhead. The table, ui(a), requires a 12 bit address; and the 
table, V2 (q,/ 3), requires 7 bits in a by 9 bits in 0 (this has not been optimized); for 
a total of 66k words. 
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Domain and Range of a and 3 

The ratio a has a. pole in the negati\*e \'elocity-domain (i.e., a red-shifted line). To 
avoid this, a test is made to determine the sign of the velocity. For positi\'e velocities: 
Q = ^ Pqy negati\’e \elocities: a = h-h ' 

The velocity is negative if and only if: > I 4 — I 2 (2.13) 



As can be seen from figure 2.3, a covers a more-or-less symmetric range about 
zero. Therefore, it can be scaled-up to provide an address for a look-up table without 
having to subtract an offset . 

The function is s\-mmetric in the velocity domain and is continuous for sym- 
metric line profdes. For asymmetric lines (and noisy lines) 3 has a pole near u = 0. 
In practice this does not cause prol)lems because, near v = 0, v is independent of 3 
(i.e., r 2 (a,,d) ^ 0 for small r). An offset may have to be subtracted from 3 when 
scaling it for use as an address. 

Finally, note that: 



l2-h ^ (/l-/3)-(/4-/2) 

(A -A) 



( 2 . 14 ) 



Hence, both a and 3 eaii be ol)tained from the differences (/i — /a), (7i — 74 ), and 

(A -A). 



The Algorithm 

It is assumed here that the difhn'ences: A\ = I\ — I 3 , Ay = I\ — I 4 , and A 3 = I 4 — I 2 
have been accumulated in image buffers. 

For all pixels: 

Calculate 3 and sign of wlocity: 

nume rator = Ai — A 3 

{ — 1 . if numerator > 0 : 

+ 1 . otherwise. 
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Calculate q: 




Calculate velocity: 



if sign > 0: 

otherwise. 



(. = sign -[11(020) + \ 2(C30. Oh /i)] 



Oi , . . . , O 4 are constants; \ i and V 2 are look-up tables. 

It should be noted that noise in the data may cause a ^'division by zero" error in 
the algorithm. 



Constructing the Calibration Tables 

The following i)rocedures are used to construct the calibration tables \\ and 1 ^ 2 - 
Parameters used for this study appear in square brackets. 

The table Vh(o) is made for a line-profile of nominal width [130 mA] as follows: 

1 . Calculate a for a representativ^e number [44] of discreet (evenly spaced) 
\'elocities ranging between some number slightly less than zero [-40 m/s] 
and some number slightly larger than the largest velocity expected from 
the data [4500 m/s]. 

2. Invert the table of n versus v to get a table of v versus a. 

3. Interpolate (cubic spline) between the endpoints to get Vi(of), a table of 
\'elocities for an evenly spaced abscissa a. 

The table made for tlie same range of velocities as is Vi; and, for a 

given range of widths [100 to 295 mA]. The minimum width should not be less the 
measurement spacing [S5 mAj. 

There are two 'Uinphysical” regions in the table I 2 - The highest values of a cannot 
occur with the highest \'alues of and the lowest values of a cannot occur with the 
lowest values of J. These portions of the table will not normally be addressed by the 
Doppler algorithm. 

V 2 is calculated as follows: 
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1. Calculate o and 3 for a represent ati\e grid of velocities [22] and widths 

[14]. 

2. Interpolate (spline) to get a much finer grid [512 in velocity by 256 in 
width]. 

3. For each “real’* V(docily v [512 of them] calculate the v^elocity correction 
C2: 

V'2 = V — \ \ (a) 

using the table \ i(o) computed abo\'e. 

1. Sort the ordeif'd triplets (cj.o, J) into a 2-dimensional table V 2 {cx,^). keep- 
ing track of how many times each bin gets filled. 

5. Fix uncoop(M'al i\'(' bins: 

a) If a l>iii accumulated more than one value of ^ 2 , normalize to get 
the “axerage" C 2 for that (o,/3), 

b) If a bin in a "physicaT' region did not get filled [there were 6 of 
them], interpolate [linearly] using adjacent bins. 

c) Fill “uiiphysicar' bins with the value of Co from the nearest ‘'physi- 
cal" bill, riiis is done because noise in a measurement may produce 
an ordered pair (o. .i) which does not lie Cjuite where it should in 
the table. 

Tables \ i and \ 2 ^ven^ made in units of m/s, and rounded to the nearest integer. 

2.2.2 Spectral Simulation 

The instrument and Doppler algorithm were simulated in order to determine the 
magnitude of statistical and systematic errors. The simulation program stored the 
line and transmission profiles in single-precision 512-element tables covering the range 
±800 mA about nominal line-center. The intensities, /j, . . . , / 4 , were calculated for a 
discrete set of x^elocities as follows: 



1. Offset the index of the line-profile table to simulate a Doppler shift. 
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2. Take the product of the transmission-profile and the offset line-profile. 

3. Integrate [sum] o\'er wav'elength to get "measured" intensity. 

The Riemann inter\ al and limits of integration were varied to verify that they 
did not affect the results. 

The error estimates are described in the next two sections. 

2,2.3 Sources of Systematic Error 

Line-Depth Variations 

The ratio a will l)e sensitive to variations in line depth if measurements H and /3 (or 
I 4 and I 2 ) do not integrate the same amount of continuum. 

Figure 2.5 shows the transmission through the Lyot and Michelson filters as a 
function of wavelength for the first and third measurements. The ratio of integrals of 
these curves is plotted as a function of Lyot filter fwhm (assuming a Gaussian Lyot- 
profile) in figure 2.6. The two integrals are equal when the Lyot fwhm- is ~ 625 mA. 
They are different by 2 % when the fwhm is 500 mA. 

This 2% difference causes a \*elocity error tywhich depends on line-depth as shown 
in figure 2.7. If the L}'ot fwluii were 400 niA, would be larger by a factor of 4. 

Lyot Filter Drift 

If the Lyot filter profile drifts in wavelength space, the four measured intensities may 
not integrate the same amount of continuum. This would give erroneous values for a 
and /? and, hence, e. To (diminate this source of error the transmission profile 
should satisfy: 




^ = 1 4 



(2.15) 



for some small range d\ over whicli the Lyot filter might drift; 1^^^ trans- 

mission of the Michelson filter at wavelength setting /: 



For a Gaussian Lyot-profile, it turns out, this condition is more or less satisfied 
for the same fwhm that minimizes due to line-depth variations. The reason can be 
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Figure 2.5: Transmission product for Lyot and Michelson filters vs. wavelength. <i 
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Figure 2.7: Systematic error vs. velocity for line-depth increments 

- 20 ,- 10 , 0 , 10 , 20 %. 



seen from figure 2.5. For transmission profile, if the Lyot profile shifts to the left 
(right), any increase (decrease) in the area under the main-lobe is compensated by a 
decrecise (increase) in the area under the largest side-lobe. Transmission profile <3 is 
unchanged to first order since it lies close to the maximum of the Lyot- profile. 
Figure 2.8 shows for S mA offsets of the Lyot filter. 

Line-Width Variations 

As stated above, if a alone is used to compute v, an error will be introduced which 
depends on the line-width. Figure 2.9 shows the error incurred if the line-width is 10, 
20, and 30% larger than the nominal width for which the calibration table Wcis made. 
If both a and /3 are used, the error is reduced to a few m/s cis shown in figure 2.10. 
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Figure 2.8; Systematic error vs. velocity for Lyot drifts; —8,0.8 mA. 




(a-only method). 
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Figure 2.10: Systematic error £„ vs. velocity for line-width increments: 0, 10, 20, .30% 
(a, method). 

Line- Asymmetries 

The effect of line-asymmetry was e.xplored by using a line-profile with a Gaussian- 
shaped bisector as shown in figure 2.11. It was found that the depth at which the 
velocity is measured depends on the velocity, but is independent of bisector shape. 
The velocity dependence is plotted for different bisector shapes in figure 2.12 (a-only 
method), and figure 2.12 (a. 3 method). Note that does depend on the precise 
shape of the bisector, but its magnitude cannot be determined until a measurement 
is made of the bisector shape at the spatial resolution of the MDI instrument. 

2.2.4 Statistical Error 

Statistical errors were found by a standard propagation of error calculation. Deriva- 
tives of the functions ei(o) and V 2 {o:. 3) were approximated by finite differences. 




depth of V(mdi) 
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Figure 2.12: Depth of measured velocity vs. velocity for bisector shape hwhm: 
—7.5, .... 7.5 mA (a-only method). 



Figure 2.11; Geometry of modeled line asymmetry. 
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Figure 2.13: Depth of measured velocity vs. velocity for bisector shape hwhm: 

—7.5, . . . , 7.5 mA (q. 3 method) 




For a sinusoidal line-profile with delta-function transmission-profiles, the statisti- 
cal error can be written down in a form like that shown in the SOI/MDl proposal. 

_2cs 1 

d 

where, s is the measurement spacing; d, the line depth; and, Iq, the number of photons 
collected at continuum intensity (half the well depth for the CCD). 

Numerically, from this calculation is 16 m/s (with d = ..58, Iq = 10®). The finite 
spectral resolution, howe\ er, tends to wash out the signal. For Gaussian transmission- 
profiles {fwhm = 100 mA). is 25 m/s. For the transmission-profiles in figure 2.5, 
<Tv is 30 m/s as shown in figure 2.14 (a-only method) and figure 2.15 (o, (3 method). 
The scatter in figure 2.15 is due to the method by which the calibration table V2(a,/9) 
is calculated. 

The sudden increase of cr^, for velocities above 3500 m/s is due to the increasing 
dependence of v on 3 together with the fact that 3 is only allotted 9 bits in the 
calibration-table address. .Allotting more bits to 3 would proportionately reduce this 
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Figure 2.14: Statistical error cr^ vs. velocity for different line-widths (a-only method). 
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Figure 2.15: Statistical error <Ty vs. velocity for different line-widths (a, /? method). 
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behavior. 



2,2.5 Discussion of the Method 

The studies indicate that velocity errors due to line-depth variations and Lyot filter 
drift can be minimized by proper design of the Lyot filter profile. In principle it 
should be possible to eliminate the errors entirely, though only a crude optimization 
has been done here. 

Line-width variations can produce velocity errors of up to a few hundred m/s 
when the velocity is calculated from a single rational function of the data (a-only 
method). The error is reduced to a few m/s by using two rational functions (a, j3 
method). The more accurate computation requires an extra image buffer and entails 
doing 2,2 times as many operations. 

A line-asymmetry will affect the measured velocity because the line-bisector is 
sampled at different depths across the solar disk. This happens because the average 
velocity is a function of disk position, and the data tend to sample more of the 
line-wings for large Doppler shifts. 

Data taken at La Palma in 19SS show that the averaging of subresolution velocity 
elements produces a “'noise'' le\’el of ~ ±20 m/s. Hence, it seems reasonable to set 
this as an upper limit for systematic errors. 

The statistical error for a four point measurement is ~ ±30 m/s with a CCD 
well-depth of 2 x 10^ photons. The error becomes significantly larger at ^ 3500 
m/s for the a, j3 method. This can be rectified by use of a larger calibration table. 
The errors for 2 pass (S point) and 4 pass (16 point) measurements are 21 m/s and 
15 m/s, respectively. 

In order to minimize the statistical error, sidelobes should be kept to a minimum. 
Hence the optimum Lyot filter-profile would have a flat peak with wings just large 
enough to ensure integration of equal amounts of continuum at each wavelength 
setting. 

This recommendation for the Lyot filter width has been used in the design of the 
MDI instrument. The choice of velocity calculation method is still pending further 
study. It is not \'et known whether line-width or line-asymmetry variations will 



CHAPTER 2. MEASUREMENT TECHNIQUES 



30 



contribute more to measurement errors. Only if line-asymmetry variations are small 
can the a. R method used. 



2.3 Magnetograms 



2.3.1 Inverse Zeeman Effect 

Spectral lines of atoms in a magnetic field are split by the Zeeman effect. For a field 
component //|| pointing toward the observer, an emission line of wavelength Aq is split 
into two components {a coiii[)onents) with opposite circular polarizations. The left 
and right circular polarizations (LCP and RCP) are displaced from Aq by — Aa and 
+A\, respectively. 

Hu 



Aa = - c/.\n (2.17) 

where g* and g” are the Lande factors, and 17' and A/" are the magnetic quantum 
numbers for the initial and final states (Bray k. Loughhead 1964). The magnetic 
quantum numbers are subject to the selection rule A/" — A/' = ±1. For the Fe I 
(6302.508 k) line. 

" (2.18) 



-—2 = 17.4 m.4/kgauss. 

OI 

For absorption lines, the LCP and RCP components will be preferentially absorbed 
at Aq — A,\ and Aq -f A,\. respectively. This is termed the inverse Zeeman effect. 



2.3.2 Longitudinal Field Measurement 

From figure 2.16 we can write down the measured LCP and RCP intensities. 



— fo + 



dl 

d\ 



Aa 



(2.19) 



Ir = ^0 — 



di 



d\ 



A > 



( 2 . 20 ) 



where ^ is assumed to be constant over the range X^j ± Aa- 
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We can now define the magnetogram signal M, 



M = 



Il-Ir 



1 



II + Ir Iq 
Combining this with equation 2.18, we have, 

//|i = cM 



dl 

dX 



A, 



( 2 . 21 ) 



( 2 . 22 ) 



where the constant of proportionality c depends on the profile of the absorption 
line and of the filter through which the intensities are measured. The constant is 
independent of continuum intensity. For the filter parameters in table 4.1, 

c % 10^ gauss. (2.23) 



2.4 Transform Methods 

The goal of the data analysis is to measure the dispersion relation and amplitude 
of waves. The natural method for plane wave decomposition in a linear medium is 
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the Fourier transform. ( 'onsider a monocluomatic plane wave p with wave vector 
^'o = (^OrN'oy) and frequency 

p{x.t) = cos(A^o • T - ^'oO (2.24) 

The Fourier transform of f> is defined as. 

Iff cos{k-Q ■ .V — dt (2-25) 

Evaluating the integral. 

~ ~ ~ ^’o) + ~h k'o)S{Lj + cc’o) (2.26) 

The two delta functions, l> ing o{>posite one-another in space, illustrate a general 

property of Fourier transforms. For real valued functions, the (complex) Fourier 
transform will have a factor of 2 redundancy in information. In fact, we only need 
to look at a half-space (containing the origin) to obtain all the information about 
the original function. It is conceptualh’ easiest to eliminate the negative frequency 
half-space, and consider only: 

— oc < k‘r < zc. — cc < ky < oc , 0 < u; < cxo 

It is important to account for relative motion of the observers field of view with 
respect to the obser\*ed medium. If the medium is moving at constant velocity v with 
respect to the observer, we have x x — vt. This produces a measured frequency 
shift u: ^ for the plane wave, where 

A^ = A- c. (2.27) 

A final and extensive!}- used property of Fourier transforms is the convolution 
theorem. It states that the Fourier transform of the convolution of two functions 
/(x) and g(x) is equal to the product of the Fourier transforms of the functions. 

nosi-r -()<!(} = X{fU)} 



(2.28) 
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2.4.1 Practical Considerations 

The discrete Fourier transform, or Fourier series, can be viewed as an extreme example 
of a Fourier transform (Bracewell 1986). The series representation Et is useful for 
discrete, evenly spaced data /„ with periodic boundaries. It is defined by, 

a = 1 = 4 E C2'29) 

n=l 

The inverse transform is defined by. 

L = (2.30) 

k = l 

For measurements made at an interval d. the maximum Fourier frequency that can 
be represented is the Ny([iiist freciuency. 

/nv, = ^ (2.31) 

The discrete transform has the ad\*antage that it can be evaluated using a very 
efficient algorithm, the Fast Fourier Transform (FFT) algorithm. The data from the 
series of Doppler and intensity images is indeed binned discretely on an evenly spaced 
grid in the spatial dimensions (by pixel) and temporal dimension (by image number). 
The condition of periodic boundaries must be enforced, however, by apodizing (see 
below). 

There are several algorithms a\*ailable for discrete Fourier and related transforms. 
For handling real valued data, it is more efficient (by a factor of 2) in computing 
time and memory usage to construct the Fourier transform from a Fast Sine-Cosine 
Transform. For doing convolutions to smooth and filter data, it is most efficient and 
conceptually simplest to use the Fast Hartley Transform. The abbreviation FFT is 
used generically hereafter to refer to the Sine-Cosine and Hartley transforms. 

The FFT algorithms used were designed for data series containing N points, where 
N is even and factorable into low primes. The spatial dimensions were 512 pixels in 
both X and ij (A'V = \y — 2‘^). The temporal dimension was 4S0 images (iV^ = 2^-3-5). 
This number was chosen because it entailed the largest ''nice” number of filtergrams 
that would fit on a single data tape. 
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2.4.2 Apodization 

I 

To approximate periodic boundaries, the data is multiplied in each dimension by a i 
function which tapers to zero at the edges. The exact form of the function must ^ 
be chosen very carefully to avoid disturbing the Fourier spectrum. The convolution 
theorem tells us that multipl\'ing the data by an apodizing function is equivalent to 
convolving the Fourier sp('ctrum with the Fourier transform of the apodizing function. 
Therefore, the apodizing function should be constructed such that its transform has 
a narrow bandwidth and \’er\’ small sidelobes. 

Several standard apodizing functions are listed in Harris (1978). The 3-term 
Blackman- Harris function was choscni. Its discrete Fourier transform has a 3 dB ‘ 
width of 1.56 bins, and largest sidelobes of —61 dB. The windowing function for N 
bins is given by 

■ do + "I cos (^'*) + (E cos (^'^) (2.32) 

tor n = — ^ + 1, . . . , y — 1 



where oq = 0.44959 



= 0.49364 



a? = 0.05677 



When filtering data <is done in chapter 7, it is not, in principle, necessary to 
apodize the data. For tlie La Palma data sets, however, the image rotates during 
the observations. Removal of this rotation from the images leaves the corners of the 
images without real information. To minimize round off errors from the corners, the 
data were masked before filtering. The mask used had a cosine taper (to zero) at the 
edges of a flat central disk: 



1 if rjriji < ^ - 32 

\ 0 if T'jYiJl ^ 

i cos - -V + 32)] otherwise 



(2.33) 



w here r. 



= \/ m- + /?*. 



f-1 



//?, n 



•2 
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2.4.3 Normalization 



The FFT’s are normalized such that the spectral power of a sine wave is ecjual to the 
square of the wave-heighl before apodization. 



“ 2 
|^(l Fn • 5 sin A^onj 

a-=-4 



= 1 



(2.34) 



The e.xact normalization is a weak function of k and is different for the real and 
imaginary parts of the transform. However, it can be approximated by a constant. 



= A\J^{W'Jn}\^ (2.35) 

(2..36) 

This approximation is valid for wavelengths that are short compared to the scale 
of variation of the apodizing function W'n^ is the case for all wavelengths and 
frequencies of interest. 

There is an additional wa\eiuimber-dependent normalization factor that corrects 
for the attenuation of spatial frec[uencies by the atmosphere and optics. This correc- 
tion is described in sections 3.1 and 5.4.1. 



wliere 



.4 = 









jV 
■ 2 



2.4.4 Statistical Error Estimates 

An error estimate can be obtained for the spectral power |Ffcp by a standard propa- 
gation of error calculation. If we assume a constant error ay on the data /„, we can 
use equation 2.29 to find I he error on the power 

(2.37) 



In three dimensions, this becomes: 



Chapter 3 



Telescope and Instrumentation 



3.1 Telescope 

The 50 cm Vacuum Solar Telescope at the Swedish Solar Observatory is shown 
schematically in figure 3.1. A photograph of the tower at Roque De Los Mucha- 
chos on the island of San Miguel De La Palma (Canary Islands, Spain) appears in 
figure 3.2. The telescope is a refractor with an / number of 45. 

The image forming lens is a 50 cm achromatic doublet. Two flat mirrors Ml and 
M2 in the altazimuth turret guide the beam down an evacuated tube. The first lens 
of the achromat forms tlu' vacuum window. At the bottom, the flat mirror M3 directs 
the beam through a vacuum window to the (Ursies) observing room. Extreme care 
was taken by the designers of the optics to ensure the telescope would perform at the 
diffraction limit (Scharmer et al. 1985). 



3.2 Postfocus Instrumentation 

3.2.1 Tunable Filter 

Lyot Filter Operation 

A Lyot filter is construct('d IVoni a sei ies of simple birefringent elements. A represen- 
tation of one such element is shown in figure 3.3. .An incident pulse of light is polarized 
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Figure 3.1: Schematic diagram of the Swedish 50 cm Vacuum Solar Telescope. From 
Scharmer et al. (19S5). 

at 45° to the optical a.\es of a birefringent crystal. Fast and slow pulses emerge from 
the crystal, and pass through an exit polarizer. The output of the element is, thus, 
a pair of identically polarized. ec[ual amplitude pulses, with a temporal separation 
proportional to the crystal thickness. The spectral response of such an element is 
given by the Fourier transform of its impulse response. The Fourier transform of a 
pair of delta functions is a cosine. 

-A Lyot filter consists of a cascade of elements, each double the length of the 
previous one (see figure 3.4). The transmission of each element as a function of 
frequency is shown in figure 3.5. The longest element has the most closely spaced 
maxima. The next longest has minima that fall on every second maximum of the 
longest element, and so forth. The net effect is to produce a transmission profile 
with a mainlobe width approximately equal to that for the longest element; and a 
free spectral range, or distance between peaks, equal to the period for the shortest 
element. 
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Figure 3.2: Swedi.sh Solar Observatory. La Palma. Spain. 
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Figure 3.3: Pulse propagation through a simple element. From Title & Rosenberg 
(1981)( 1 




Figure 3.4: Lyot filter with four elements. From Title Rosenberg (1981). 
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Figure 3.5: Transmission vs. frequency for the components of a Lyot filter (A, B, 
C, D) and the transmission through successive stages of the cascade (E, F, G). From 
Title Rosenberg (19S1). 
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Figure 3.6: Birefringent tuning element. From Title & Rosenberg (1981). 



Tuning a Birefringent Element 

The bandpass of a birefringent filter can be tuned over the full free spectral range 
by individually tuning the component elements. The best practical method for doing 
this is illustrated in figure 3.6 (Title Rosenberg 1981). As in figure 3.3, a fast 
and slow pulse emerge from the crystal parallel to the feist and slow axes. The delay 
between the pulses is in\'ersely proportional to wavelength. Reconstituting the two 
polarizations gives elliptically polarized light. The axes of the ellipse are horizontal 
and vertical, and the ratio of the axes, the ellipticity, is a function of wavelength. 

A quarter wave plate transforms the polarization from elliptical to linear, and the 
orientation of the resulting polarization vector is dependent on ellipticity. Finally, 
a rotatable half wave plate and fixed exit polarizer are used to select the linear 
polarization, and hence, the wavelength. 
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The Solar Optical Universal Polarimeter 

The Solar Optical Universal Polarimeter (SOUP), developed by LPARL, employs a 
nine element tunable birefringent filter. It is similar in design to a Lyot filter, but each 
birefringent element contains a partial polarizer to apodize its transfer function. The 
result is a filter with smaller sidelobes than a Lyot filter. The SOUP instrument was 
used on the Spacelab 2 shuttle flight, and the filter has become the evaluation model 
for the Orbiting Solar Laboratory — Coordinated Instrument Package (OSL-CIP). 
Software-selectable blocking filters, each 6-8 A wide, are centered on solar absorption 
lines. Waveplates can also be selected to measure left and right circular polarizations 
as well as linear polarizations spaced e\’ery 45°. 

At a wavelength of 5S96 .\. the free spectral range of the tunable filter is 14 A, 
and the/u’//m of the passband is 62 mA. A *'wide mode" is also available by rotating a 
half wave plate within th<' longest element. This effectively shortens the element by a 
factor of 3, and increases the bandwidth to 103 mA. The stepping motors that rotate 
the tuning half wa\'e plates permit timing in 3.14 mA steps. Tuning is repeatable to 
within half a step. 

Temperature gradients and imperfections in the filter components produce small, 
but measurable, variations in filter characteristics over the field of view. The center 
of the passband varies b\’ up to 10 mA and the fwhm, by up to 5 mA over the field of 
view. In addition, the tobil amount of light transmitted (Integrated over wavelength) 
is modulated as the filter is tuned. (See figure 3.8.) This modulation accounts for 
most of the error in the \elocity and intensity measurements. 

3.2.2 CCD Camera 

The CCD (charge coupled de\4ce) camera uses a Texas Instruments 1024 x 1024 pixel, 
virtual phase detector. Il is also an e\^aluation model for OSL-CIP. Subareas of the 
full detector can be read out to reduce the time required to write data to tape. In 
addition, 2 x 2 groups of |)ixel.s can be summed to decrease exposure time by a factor 
of 4. This 2 X 2 summing mode was used for the majority of observations listed in 
chapter 4. 
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Figure 3.7: Transmission profile of the SOUP tunable filter vs. wavelength. Trans- 
mission was measured using a spectrometer, with the filter tuned to 6328 A. 

The detector geometry consists of a rectangular grid of pixels. Each pixel is 
efficient at collecting light up to the boundary with its adjacent pixels. Thus, the full 
area of the CCD array is sensitive to light. 

The charge accumulated on a pixel is proportional to the integrated light falling , 
on that pixel. The signal also integrates a current present in the absence of light, i 
the “dark current." The dark current depends on device temperature. To keep this ' 
effect constant, and to a minimum, the CCD is kept in a vessel at low pressure, and 
is cooled thermoelectrically to —20° C. 

3.2.3 Image Stabilization System 

Compensation is made for small image displacements, caused by seeing and inaccurate j 
telescope tracking, by an active tilt mirror in the beam path. The mirror is moved 
in accordance with signals from a sunspot tracker which senses image motion. The 
active mirror is mounted on three piezoelectric actuators which can tilt the mirror on i 
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Figure 3.8: Relative Iran.'iinission of the SOUP tunable filter vs. tuned wavelength. 
Transmission was measuiecl using a continuum illumination source. 

two axes. The system has a bandwidth of 65 Hz. 

The spot tracker is an array of four photodiodes (a quad-cell). A sunspot image 
centered on the quad-cell will produce equal signals from each diode. If the image 
moves toward one of the four quadrants, however, that quadrant will experience less 
illumination and produce a smaller signal. The differences of signals from the four 
quadrants are used to determine compensatory signals for the active mirror. The 
integrated compensatory signals are used as a feedback to the telescope tracking 
system. 



3.3 Optical Layout 

The optical bench was set up in the Ursies room as shown in the photograph in 
figure 3.9. The schematic layout of components is shown in figure 3.10. The beam 
from the telescope strikes the active mirror, and is deflected through a field stop 
slightly larger than the active area of the CCD. The heat stop, used to filter out 



CHAPTER 3. TELESCOPE AXD IXSTR [MENTATION 



45 



infrared wavelengths, was removed from the light path because it produced fringing 
patterns. 

Wavelengths of 5000 7000 A are reflected by the dichroic beam splitter to the^ 
tunable filter path. Infrared wavelengths are not reflected. Other wavelengths are^ 
reimaged through the Ca 11 Kdine filter and the 4308 A bandpass filter to the K-line 
and '‘white light" cameras, respective!}'. The white light images were used to monitor 
the seeing conditions and the focus in real-time. The K-line camera was not used for 
these observations. 

After reflection from the dichroic beam splitter, the beam is collimated by a lens 
on a motor dri\*en stage. The magnification of the system is chosen by substituting, 
lenses of different focal length. The stage is moved automatically to compensate for 
the chromatic dispersion of the telescope, and the whole assembly, moved manually to 
adjust the focus. The beam passes through an aperture stop and is reimaged through | 
a set of selectable wa\'eplates and blocking filters, and the tunable filter. For solar- 
observations, the beam passes through a sector shutter to the CCD camera. To check 
the position of a spectral line, the flip mirror was inserted into the beam to deflect it 
to the photometer. 

Part of the beam is diverted, before the polarizer wheel, to the spot tracker. There 
is a neutral density filter and adjustable zoom lens system in this path to optimize) 
the intensity and magnification of the image for the spot tracker. 

The motor driven stage, polarizer wheel, blocking filter wheel, and tunable filter f 
were controlled by a (PC compatible) COMPAQ computer. The shutter and CCD^ 
camera were controlled b\ a \'AX .Station 3200. The system was free-running, with no 
timebase to control the data acquisition. The timing of the images, therefore, varied j 
as shown in figure 3.11. The time mark given to each image came from a crystal! 
controlled clock, and is accurate to 1 ms. 
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Figure 4.6: Optical bench setup. 
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OPTICAL LAYOUT AT SWEDISH SOLAR OBSERVATORY 
LA PALMA, SUMMER 1990 
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Figure 3.10: 
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Figure 3.11; Timing offset for images vs. time from start of run. The offset is the 
difference between the image time mark and a straight line fit to the time marks. 

3.4 Modulation Transfer Function of the System 




The modulation transfer function (MTF) of an optical system describes the effect of 
aberrations and diffraction on an image formed through the system. If Ig{x) is the 
image intensity distribution e.xpected from aberration-free geometrical optics, then 
the real image intensity distribution Ir{x) is given by, 

IAt) = J^-‘{MTF(1') . (3.1) 

The attenuation of the spectral power is given by the square of the MTF 

functions. 

The components which contributed significantly to the total MTF are the entrance 
pupil, the detector pi.xel shape, and the earth’s atmosphere. The total MTF is the 
product of the MTF’s of these individual components. The atmospheric MTF will be 
treated in section 5.4.1. 
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3.4.1 Entrance Pupil 



The entrance pupil is defined by the diameter d (4S cm) of the heat shield diaphragm 
in front of the objecti\’e. Idie MTF of the pupil is giv-en by the autocorrelation of the 
pupil function (a disk). For wavenumbers, 0 < A’ < A\Vyq^ 



MTF(A0 = 



cos ^ K — K \/l — 



.•2 



(3.2) 



f rf pixels \ k 
2(1 \ field of viewy A\\yq 



3.4.2 CCD Camera 

The effect of summing liglit o\*er tlie area of a CCD pixel is equivalent to convolving 
the image distribution with the pixel shape, and sampling the new distribution at 
points centered on each pixel. Thus, the MTF of the detector is the modulus of the 
Fourier transform of the pix(4 function (a square). 



M fF(A‘) = sincK.j: sincK.^ 



(3.3) 



1 ’ 1 

where k = 



2 fess; 



and siiic.r = 



yq 

Sm 7TX 
7TX 



This can be approximated by a cylindrically symmetric function with no more than 
a 2% error. 

MTF(A-) = sincK (3.4) 



The square of the MTF functions are plotted in figure 3.12 for the entrance pupil and 
CCD detector. 



MTF 
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Figure 3.12; Squared MTF of telescope and CCD detector 



Chapter 4 
Observations 



The original topic of this dissertation was “sunspot seismology,'’ that is, the study 
waves propagating into and out of sunspots. That topic was set aside when waves were 
discovered at substantial!}' higher wavenumber and frequency than had previously 
been observed. The new topic of anal\*sis became the characterization of properties 
of these wa\’es. Man}' of the observing runs are, therefore, extraneous to the analysis 
described here. They ha\ e becui included in the list of data in this chapter, but no 
calibration or reduction has Ix^en pen formed on them. 



4.1 Observing Procedure 

The observing runs consist of a re|)eating sequence of filtergrams with a temporal 
spacing of about 3.5 s between filtergrams. The wavelength, filter width, polariza- 
tion, and relative exposuie time for each filtergram in the sequence is defined in the 
sequence file residing on the (’OMPAQ computer that runs the filter. 

Calibration data was taken before and after observing runs. The order of calibra- 
tion and observations was as follows: 

1. Dark Currents: 4 measurements of dark currents for each exposure setting. 
Done at the Ix'gi lining and end of each day. 
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2. Line Centering: Delinilioii of the nominal line center from which wave- 
length offsets are measured. The photometer detector was used to inte- 
grate light o\’er the field of view with the target region centered in the 
field. 

3. Slide Projector Flats: t repetitions of flat field measurements using slide 
projector sourc e. 

4. Observing Run: One or more tapes-worth of observations of target region. 

5. Solar Flats: 4 repetitions of flat field measurements using defocused sun- 
light source. 

6. Slide Proj(‘(4oi Fhits: R('[)eat of step 3. 

The line centering st(‘p is performed by scanning the filter through a range of 
wavelengths (from short to long), and plotting intensity versus wavelength. This step 
was not performed correc t Iv due to a design flaw in the photometer interface with 
the COMP.AQ computer. The computer did not wait for the full signal integration 
time before reading the light intensity. In the resulting line profile scan, the intensity 
points lagged behind the wav'elength points. This caused the line to appear shifted 
to longer wavelength. 4 he line was, therefore, "centered" to an erroneously short 
wavelength. This offset of — 10 m.A for the Na Di line was discovered by making a 
line scan using the CCD ( am(u*a instead of the photometer (see figure 4.1). 

4.2 Calibration Method 

The purpose of calibration is to produce images with a signal to illumination ratio 
that is independent of both the pixel location and the wav’elength setting of the 
tunable filter. The procedure described here is more complicated than the standard 
one used by the Lockheed group. I'his is necessary to provide an accurate calibration 
for dopplergrams. The flat fi(4ding procedure calibrates the whole system, optics and 
detector together. 

The CC’D is a v’erv’ liii(‘ar ch^’ice. provided the light intensity stays below saturation 
levels. Response of a pix(4 can l)e characterized by a gain and an offset. The offset 
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Figure 4.1: Line scan of N'a D] Line using the CCD camera. Signal, integrated over 
the field of view, vs. wavelength offset. 

is determined by the amount of charge that accumulates on the CCD in the absence 
of illumination (the ‘■‘dark current"). It must be measured for each exposure time. 
Each pixel has a characteristic gain and dark current. The tunable filter is also 
responsible for some variation of the signal strength as a function of wavelength A 
and spatial position x. The net response of the system is to produce filtergrams F,(x) 
at wavelengths A,, given a source intensity distribution 5(x, A). 

Fi(x) ■= J S{x. X) A,{x, X) d\ + D{x) (4-1) 

where i4j(x. A) is the product of filter transmission and CCD gain at filter setting A,, 
and D{x) is the CCD dark current. 

Because the filter bandpass is narrow compared with spectral features in the 
source, we can approximate: 

J S{.T, X) A,{x, X) dX = S{x, Xi) jAi{x,X)dX. (4.2) 




Thus 



F,(.r) = .^fx. A,) B,{x) + D{x) 



(4.3) 
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wliere B^(x) — J A^{x, A) dX. (4.4) 

In order to clcterniiiie S’, the quantity of interest, \ye must first measure and D. , 
The dark current D is easily measured by making filtergrams with no illumination 
source. For S' = 0, we ha\e D = 

If we had an illumination source that were both flat in spectrum and flat across' 
the field of view, we could measure directh'. The sources available, however, were: 

1. Defocused Sunlight S'.s(A): Flat across field, but not spectrally flat. 

2. Slide Projector Laiu|> S’/ST): Flat in spectrum, but not flat across field. 



The technique is to make uiie filtergram Ff at wavelength A^ in the continuum ’ 
(away from the absorption line), using defocused sunlight. And, to make filtergrams 
at Ae and at A^ using the slide projector lamp. This gives us the following ' 
three measurements: 



F^^{x) — Ss{Xc) Bc{x) + Ds{x) 


(4.5) 


q-(.F) = SU-r) IM-r) + DU-r) 


(4.6) 


/•q(,r) = B,{x) + Dl{x) 


(4.7) 



Giv'en the dark currents and 72^, we can solve for the gain correction Qi at A^ 
defined as Qi(x) = Bi{x)~\ 



^ If,'- - Dl\ 

[t't - DLllff - Ds] 



(4.S) 



There is some arbitrary cA'erall imrmalization which depends on Ss{^c)* Iri practice 
the spatial average of at the A^^ filter position was defined as unity, and the other 
gain corrections, normalized accordingly. 

With the gain corrections and dark currents in hand, we can determine arbitrary 
source distributions from filtergrams. 



5(.F.A,) = c;,(.r) [C(.?) - D(.r)] 



(4.9) 
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4.3 Solar Observations 

The sequence of filtergrains for the sunspot seismology study was named SEISNAFE. 
It was designed to determine the following quantities: 

1. Doppler shift in mid-chromosphere (in sunspot umbrae). 

[II I (Balmer) (\ — line core] 

2. Doppler shift in photosphere. 

[Na D] — wings] 

3. Doppler shift in low cliromosphere. 

[Na Di — core] 

1. Doppler shift in photosphere. 

[Fe 1 (5576 A)] 

5. Continuum int{msily in green and red. 

[5575.A, 5S9 1 A] 

6. Photospheric magnetic field (line of sight component). 

[Fe I (6302.A)] 

The sequence repeated e\ery 17 s. Filter parameters appear in table 4.1. The corre- 
sponding calibration seqiuMices were N AFEFLAT for solar flats, and SEISSP for slide 
projector flats. The obser\*iug runs of appreciable length are listed in table 4.2. 

The sequence of filtergrams used for the high-frequency oscillation study was 
named NAFAST. It was designed to make dopplergrams in the low chromosphere 
and photophere, using filtergrams in the core and wings of Na Di, respectively. All 
other measurements wen' eliminated to increase the repetition rate to once every 
17.7 s. The sequence parameters appear in table 4.3 and the corresponding calibra- 
tion sequences were N.Al'EFL.AT for solar flats, and NAFASTSP for slide projector 
flats. 

The NAF.AST sequence was run only once on August 2, 1990 beginning 13:51 
UT. The target region was an active region (AR61S1) containing a small, decaying 
sunspot. This is shown on a full-disk magnetogram in figure 4.2. The run lasted 
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Table 4.1: Filter 


Parameters for SEISNAFE Sequence. 


Line (A) 


Oirs(‘t 

(mA) 


Filter Mode 
(width) 


Polarization 
(0 — linear) 


Purpose 


H I o 6562,8 


-900 

0 

100 


wide 


0 


dopplergram 


Na Di 5S!)5.9 


-10 

0 

-10 


narrow 


0 


dopplergram 


Na Di 5895. 9 


- 1 75 
-125 
125 
175 


narrow 


0 


dopplergram 


\a Di 5895.9 


-2000 


narrow 


0 


continuum 


Fe I 5576.1 


-90 

-30 

0 

30 

90 


wide 


0 


dopplergram 


Fe I 5576.1 


-750 


wide 


0 


continuum 


Fe I 6302.5 


-60 

-60 


narrow 


LCP 

RCP 


magnetogram 
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Table 4.2: La Palma Observing Runs using SEISNAFE Sequence. 



Date 

(1990) 


Time 

(UT) 


Duration 

(min) 


Seeing 


FOV 


Target Region 


28 Jul 


7:47 


113 


goocl/v.goocl 


75" 


AR6168 (medium sunspot) 




10:26 


112 


lair/var. 


150" 


AR6168 


29 Jul 


8:39 


156 


fair/ var. 


150" 


AR6174 (medium sunspot) 




16:53 


126 


goocl/va r. 


150" 




30 Jul 


7:50 


281 


variable 


224" 






13:12 


293 


va riable 


224" 




31 Jul 


8:34 


237 


goocl/var. 


224" 


flare in morning 




13:09 


347 


variable 


224" 




1 Aug 


8:08 


256 


good /var. 


224" 


AR61S1 (small sunspot) 


1 Aug 


13:35 


310 


good /fair 


224" 


AR6179 (big sunspot) 


2 Aug 


7:46 


151 


good /lair 


150" 


— overexp’d for umbra! osc. 


2 Aug 


10:52 


101 


fair/ var. 


224" 


AR61S6 (big sunspot on limb) 


5 Aug 


9:20 


190 


fair/ var. 


224" 


quiet region near AR6192 



Table 4.3 


: Filter Parameters for N.AF.AST Seq 


uence. 


Line (.A) 


Olfset 

(mA) 


Filter Mode 
(width) 


Polarization 
(0 == linear) 


Purpose 


Na D, 5895.9 


-40 

0 

40 


narrow 


0 


dopplergram 


Na Di 5895.9 


- 1 75 
-125 
1 25 
1 75 


narrow 


0 


dopplergram 


Na Di 5895.9 


-2000 


narrow 


0 


continuum 
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Figure 4/2: luill-disk niagiu‘logr<un from the Xatioual Solar Observatory at Kitt Peak. 
The field of \ ie\v coiit ainiiig AH0181 is outlined. 

143 mill with good/\'ariable that deteriorated in the last 30 min of the run. To 

obtain more information al)onl the acti\'e legion, the SEISXAFE sequence was run | 
before and after tlu‘ X.M* AS 1 runs. A continuum image is shown in figure 4.3, a high ^ 
resolution magm'logram. in hguK* 4.4. and. an offl^and H a image, in figure 4.5. 
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Figure 4.4: Magnetograin of AR()181. Network field is visible around the sunspot. 
The small tick marks are 5". 
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Figure 4.5: II a image of .AR6181 (400 m.A in the blue wing). The small tick marks 



are 5". 






Chapter 5 

Data Reduction 



The raw data used for this dissertation reside on a (full) 2.2 Gigabyte Exabyte tape. In 
producing a power spectrum, the \’olume of data is reduced by a factor of ~ 3 x 10^. 
There is more to the reduction procedure, however, than extracting a manageable 
amount of information from a cuml:)ersonie data set. The quality of the data depends 
on careful calibration to account for instrument response, image rotation and drift 
caused by telescope motion, and, distortion caused by the atmosphere. The absolute 
spatial scale of the images must also be determined accurately. This chapter describes 
the details of and cross-checks for the data reduction procedures. 

The reduction steps are listed below in the order they were performed: 

1. Flat Fielding: Calibrate for variation in filter transmission and CCD gain, 
and subtract CCD dark current. 

2. CCD Flaw Correction: Interpolate over nonfunctional pixels. 

3. Correction for Atmospheric Blurring: Normalization of measured MTF. 

4. Rigid Rotation: Compensate for image rotation caused by telescope track- 
ing motions during obser\'ations. 

5. Correction for .Atmospheric ‘^Shimmer.^^ (Destretching) 

6. Calculation of Dopplergrams and Line Center Intensity. 

7. Rigid Alignment: C’ompensate for image drift caused by telescope tracking 
errors. 
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S. Fourier Transform of Dopplergrams and Intensities: Described in 
section 2.4. 

9. Angular Scale Determination: Convert units from pixels to arcseconds. 

10. Distance Scale Determination: Measurement of foreshortening and direc- 
tion toward di^k center. Described in section 6.1.1. 

11. Measurement of Oscillation Properties: Described in chapters 6 and 7. 

5.1 Flat Fielding 

The flat fielding scheme is explained in section 4.2. We can use equation 4.9 to 
calibrate filtergrams, given dark currents D and gain corrections Qi, The dark current 
is an average of four filtergrams made with a light stop placed in front of the filter. 
The gain correction is constructed using equation 4.S as follows. 

The three flat field measurements , and F^ are each made from an average 

of four filtergrams. Because the filter stepping motors did not always find the correct 
positions, it was sometimes necessary to throw out one of the four measurements 
and use only three to make the a\'erage. A filtergram was thrown out if the spatial 
average of the central 64 x 64 pixels of the image differed by more than 5(J from the 
spatial averages of subareas of the other three images, where a is defined by inefficient 
statistics, 

< 7 , = ~ Allin (5 1) 

A I 

The gain correction is the product of (F/ — Ds)~^ and {F^ — Di)/{Ff" — D^), 
The latter is the ratio of the transmittance in continuum to transmittance an A,*. 
This ratio should be slowly varying over the field of view. The dark currents are 
not constant with time, however, and tlie very long exposure times required for slide 
projector flats (2-3 s) caused accumulation of large charges at some small groups of 
pixels. It was. thus, necessary to '‘deflaw'' F^ and F^^ (see next section), and smooth 
the transmittance ratio using a Gaussian point-spread function of 18 pixel fwhm. 

The procedure described here produced (flat fielded) images with signal to illu- 
mination ratios that were constant to within about 3% across the field of view, and, 
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more importantly, to within about 5% across the wavelength range of an absorption 
line. These numbers were measured from filtergrams of defocused sunlight (the flat 
field illumination source) at continuum wavelengths. 



5.2 CCD Flaw Correction 

Of the 2.5 X 10’^ pixels (in 2 x 2 summing mode), 4513 were not functional. An 
additional 491 pixels had exceptionally high dark currents for exposure times greater 
than about 500 ms. To minimize problems caused by the missing data, values were 
filled in at these pixel locations by linear interpolation from adjacent pixels. This 
procedure is termed "deflawing*’ . 

A good map of nonfunctional pixel locations was available for exposure times 
below 500 ms. For the longer exposure times, a map was constructed using the dark 
current runs. To do this, a cubic trend was fitted to each row of a filtergram. Pixels 
which differed from the trend hy more than 7% were flagged as “bad”, and their 
locations, entered in the map. 

5.3 Image Rotation 

The altitude-azimuth a.\'(\s of minors .Ml and M2 cause the solar image to rotate 
during the day as the telescope tracks the sun. The rotational direction and rate are 
functions of date and time of day. and can be calculated from geometrical consider- 
ations. Standard software e.xists to determine the angle of rotation and regrid the 
images by bicubic interpolation. This removes the apparent rotation in a sequence of 
images, but some data is lost from the corners of the images. 



5.4 Atmospheric Distortion 

Turbulence in the earth's atmosphere will modify wavefronts of light before they reach 
the telescope. Any plane wave passing through the atmosphere will be distorted. 
The resulting wavefront, however, can be approximated by a set of planar patches. 
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The isoplanatic angle describes the angular extent over which image distortions are 
correlated. For excellent seeing conditions, it can be over 100" (von der Liihe 1988). 
For poor seeing, it can b(' less than I". 

A change in angle of incidence of a planar wavefront produces an image displace- 
ment. When dilTerent parts of an image are displaced in different directions, the 
image appears distorted. The distortions change with the atmospheric turbulence. 
This can occur ov'er a range of timescales, from 1 to 1000 ms, with most power in the 
10-100 ms range. 



5 4.1 MTF Normalization 

The distortion of a wa\'efront from a gi\’en object point differs over the telescope 
aperture. Moreover, local displacements of up to an arcsecond cause the image to 
move around in times of order 10 ms. For exposures lasting tens to hundreds of 
milliseconds, these changing distortions will be superimposed to produce a blurred 
filtergrain. The effect is equivalent to convolving the image with some unknown 
point spread function. The Fourier transform of the point spread function is called 
the MTF. It is assumed lun^e that the point spread function is the same everywhere 
in the field of \dew of a gi\'en filtergram. 

The atmospheric MTf^ can be estimated from the spatial frequency spectrum of 
a filtergram 5(/). A contour plot of spatial frequency spectra is shown in figure 5.1 
for the —40 mA filtergram set. The sharpest filtergrams in the set are assumed to 
be diffraction limited (atmospheric MTF = 1). Since the structure of solar images 
should not change appreciabl}' over the observing run, the spatial spectra of diffraction 
limited images should be a constant Any difference between the spectrum 

of a given filtergram and .Smax nmst, therefore, be due to the atmospheric MTF. 

.S'(/) = .\[TF.,„,(/).5n,ax(/) (5.2) 

In practice the log-spectra were compared with the mean log-spectrum log 5. 
Log-spectra of individual liltergrams were smoothed by convolution with a cosine-bell, 
and the difference between log S and the smoothed spectrum was used as a correction 



Log(RMS amplitude) 
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Figure 5.1; Spatial frequeno’ log-spectra of Na Di —40 mA filtergrams. The solid 
central line is the mean log-power; the dashed lines are Icr contours; the outer dotted 
lines are the maximum and minimum power at a given frequency of any image. 
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factor to the filtergram fr('queiicies. Thus, the *\MTF normalization” procedure forced 
the spatial spectrum of each filtergram to look like the average spectrum. 

In principle the indi\udual filtergrams could have been normalized to 5max- This 
would have involved very large correction factors for the poorer filtergrams, and would 
have boosted the noise in these filtergrams considerably. Noise in filtergrams produces 
noise in dopplergrams and line center intensities in a nonlinear manner. So, instead, 
the difference between the mean spectrum and 5max retained as a correction 
factor to be applied to the oscillation power. (See chapter 6.) As a final caveat, it 
should be noted that even the best images, assumed above to be diffraction limited, 
actually include some atmospheric blurring. Hence, the correction factor will provide 
an underestimate for powcu’ al the highest wavenumbers. 

The MTF normalization procedure was checked by comparing filtergrams taken 
in poor seeing with filtergrams taken at almost the same time in good seeing. After 
normalization, both filtergrams showed the same solar structures. 

5,4 ,2 Dest retching 

The displacements cause<l b\' atmospheric turbulence vary over the field of view. 
They are independent of one another over angular scales greater than the isoplanatic 
angle. Relati\'e displacements between successive images can be measured by local 
correlation tracking (No\ ember 19S6). This procedure finds the displacement of a 
small subarea of one image with respect to the other from the cross-correlation of the 
two. .A displacement \*eclor field is measured using a patchwork of subareas of the 
images. .A relative disi^lacement grid of 64 x 64 vectors was made, correlating each 
filtergram with the preceding one. This was done using standard software. 

Local displacements are present not onl}’ because of seeing, but also because of 
wave motion and horizontal flows on the sun. We wish to remove the seeing displace- 
ments, but keep the real solar motions. It is possible to distinguish the two on the 
basis of timescale of the motion. Horizontal flows have characteristic timescales of 
hours, oscillations have timescales of minutes, and seeing has a timescale of tens of 
milliseconds. For filtergrams taken seconds apart, seeing motions will appear random, 
while solar motions will l)e coherent. 
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The seeing displacement grid was determined from the relative displacement grid 
as follows. For each v^ector component at each grid location, the relative displacements 
of preceding filtergrams were summed to produce a total net displacement (from the 
first frame). The total displacement was smoothed by convolution with a Gaussian 
of 90 s (5 bin) fwhm, and this smooth trend, subtracted to produce the seeing grid. 

It is important to ha\’c high temporal resolution in the data. Correlation tracking 
can lock onto the crests of waves. So, if the impulse response function, used to smooth 
the data, is wider than 120 s, the oscillations will be removed from the data along 
with the seeing. The standard Lockheed method for finding the seeing grid is designed 
for horizontal flow studies. It uses boxcar averages of 20 or more bins. This would 
be inappropriate for seismology studies. 

With the seeing displacements measured, the images can be ''destretched” to 
their undistorted states. \ \\e algorithm to do this treats images as rubber sheets, and 
stretches the sheets (by regridding) based on the local displacement vectors. This part 
of the process also used standard software. The set of filtergrams for each wavelength 
was destretched independently of the other wavelength sets. Destretched filtergrams 
taken at the same time, but different wavelengths, were then visually compared to 
ensure their features lined up. 



5.5 Dopplergrams and Intensities 

Velocity in m/s and line center intensity in arbitrary units were calculated at each 
pixel using equations 2.2 and 2.3. Because of the line centering problem described in 
section 4.1, the line center velocity and intensity were calculated using the —125, —40, 
and +40 mA filtergrams. The actual offsets of these measurements from the average 
line center was about —85.0, and +80 mA. Each dopplergram and intensity image 
was given a time mark a\eraged from the time marks of the three filtergrams. 

Filter tuning errors and seeing, which deteriorated toward the end of the run, 
produced some poor quality filtergrams in the data set. This in turned caused some 
nonsensical velocity and intensity determinations. Pixel velocities that were greater 
in magnitude than 10 km/s and intensities that were a factor of 10 greater than 
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average were replaced with the mean velocity/intensity for that frame. 

5.6 Image Drift 

The rotational regriddiiijz, assumes the images were rotating about a central point. 
In fact they rotate about the point centered in the spot tracker. This inconsistency 
produces a translational motion of the solar image. The point (on the sun) centered 
in the spot tracker, moreox'or, changes with time. This is caused by variable weighting 
of parts of the active region as the region evolves and as seeing conditions change. 

From equation 2.27. it is evident that to keep oscillation frequencies from shifting 
by more than 60 //Hz (O.o rr(‘([uency bins) at a wavenumber of 3 Mm“^, the image 
drift must he less than 120 m/s. This is equivalent to 1.4 pixel/hr. Image drift is 
measured from the cross-correlation of successive images. (It is defined as the distance 
of the maximum of tlu' cross-correlation from the origin.) 

The image drift is traced out o\'er time in figure 5.2. It was measured using the 
central 64 x 64 pixels of (wer\' tenth —125 inA filtergram. As can be seen from the 
figure, drift can be measured to significantly better than 1 pixel. The measured drift 
was smoothed by convoh ing the time series of each component of the vector with a 
Gaussian of 6 h\n fwhnK The trend was then interpolated to the dopplergram and 
intensity time marks and remoN'ed from these images by regridding. 

5.7 Angular Scale Determination 

For full disk images of tlie sun. the image scale can be determined by measuring 
the radius of the solar disk. For high resolution images, however, the full disk is 
not visible, and indirect methods must be used. The scale of the La Palma images 
(O.Sx lens) was determined l)y compaiTson of magnetograms from La Palma and the 
Kitt Peak National Obser\'atory (figures 4.4 and 4.2). The Kitt Peak magnetogram 
was made at 15:17 UT on .August 2, 1990, and has a scale of P'/pixel to within 1% 
(Harvey 1991). The La Palma magnetogram was interpolated from magnetograms 
made at 13:35 FT and 16:27 FT. 
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Figure 5.2: Image drift during observing run. Image motion determined by cross- 
correlation is traced out for the observing run. Points represent filtergrams taken 177 
s apart; solid line shows low-pass filter of points; circles represent 18 min time marks 
(t = 0 is at upper right). 
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The interpolation method accounted for evolution of the activ^e region by measur- 
ing the motion of the nicj^iielic field structures using local correlation tracking. The 
two magnetograms were then stretched by the time- weighted local displacement grid. 
Finally, the time- weighted signals of the stretched magnetograms were added. 

To compare magnetograms, 43 points on the La Palma magnetogram were flagged 
manually by cursor on a PIXAR computer. Corresponding points were flagged on the 
Kitt Peak magnetogram. Distances between a point and the centroid of the 43 points 
were then compared for the two magnetograms. The field of view for the La Palma 
images was, thus, determined lo be (224.39 ± 0.70) Kitt Peak pixels. Assuming the 
Kitt Peak scale measurement. 



La Palma Field of \*iew = 224^' ±2'' (square) 



(5.3) 



Chapter 6 

Measurements from Power 
Spectra 



The sets of dopplergrams and intensities were apodized in the space and time di- 
rections with Blackman- 1 [arris functions (equation 2.32). Three-dimensional FFT’s 
were calculated for the full field of view, as well as for the three 79" x 79" subfields 
shown in figure 6.1. The lemporal and spatial Nyquist frequencies were 2S.3 mHz and 
9.9 Mm”\ respectively. The frequency bin size was 117 /HIz, and the wavenumber 
bin size, 0.039 Mm“^ for the full field, and, 0.11 Mm“^ for the subfields. 

The squared modulus of each Fourier transform was calculated to produce a 3-d 
power spectrum P{k\iy). The dispersion relations given by equations 1.7 and 1.9 
are symmetric about the k = 0 axis. They define "trumpet" shaped surfaces (Hill 
19SS) of approximate form // oc (See figure 6.2.) This symmetry allows us 

to integrate the power over the azimuthal coordinate, in order to increase the signal 
to noise ratio. The result ant, collapsed power spectrum P{k,i/) will show ridges of 
power at frequencies given by the dispersion relations for the p-modes and /-modes. 

Properties of the oscillations, such as the dispersion relation and mode power, were 
measured from the collapsed power spectra. This was accomplished by fitting one- 
dimensional slices of the spectra. Fits were then compared with theoretical predictions 
to evaluate models of the solar atmosphere. Magnetic versus nonmagnetic regions 
were compared to search for differences in wave properties between the two regions. 
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Figure 6.1: .Magnetogram of .\R61S1. The three outlined subfields were used sepa- 
rately to measure wave properties in network field and quiet regions. The small tick 
marks are o". 
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Figure 6.2: Surfaces of power in space. 
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6.1 Collapsing 3-d Power Spectra 

The azimuthal symmetn' of equations 1.7 and 1.9 is realized only in solar coordinates, 
comoving with the gas. In practice, one or both of these conditions may not be met. 
The frequency shift descril:>ed by equation 2.27 occurs, with a cosine dependence on 
azimuth 0, if there is a relative (horizontal) motion between the observer’s coordinates 
and the gas. The rigid alignment procedure described in section 5.6 keeps the observa- 
tional coordinates fixed with respect to the sunspot. It has been shown, however, that 
there is relative motion b('tween magnetic flux and nonmagnetized gas (Snodgrass & 
Ulrich 1990). Hill (19SS) finds flows of 100 in/s below the photosphere. 

\dsual inspection of the power spectra did not reveal evidence of flows. Thus, 
if present, flows produce fre([uency shifts of less than one frequency bin (117 /l/Hz). 
If plane waves propagatc'd isotropically, a sinusoidal frequency shift with 0 would 
serve to widen the ridges in the collapsed power spectra. In the measured spectra, 
however, power depends mildly on angle 0. Thus, a horizontal flow would produce a 
slight asymmetry in addition to broadening a ridge. 

The second condition affecting the symmetry of the trumpets is geometrical. For 
observations awa\’ from disk-center, the projection of solar coordinates onto obser- 
vational coordinates cause's foreshortening. Waves propagating toward or away from 
disk-center appear with shorter wa\'elength (higher wavenumber) than waves prop- 
agating in the perj)endicidar direction. This foreshortening distorts the cylindrical 
symmetry of the trumpets to give them an elliptical rather than circular cross section. 
Procedures to measure and compensate for this effect are described below. 

6.1.1 Distance Scale Determination 

The distance scale of the images or, ecpiivalently, the wavenumber scale of the power 
spectra follows from trigonometry, once the magnitude and direction of foreshortening 
have been determined. This measurement was made by fitting an ellipse to the cross 
section of a trumpet in a d-d power spectrum. 

A constant freciuency slice, in the 5-miiiute band, of the full-field Doppler spectrum 
is shown in figure 6.3. The slice has been remapped to cylindrical polar coordinates. 
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Figure 6.3: Constant t'iec|iiency slice of the 3-d power spectrum. A slice at 3.99 mHz 
has been remapped to cylindrical polar coordinates. Log{poicer) is sliown as a function 
of wa\’enumber and azimuth. 

So, the cross section of a trumpet ap[)ears at roughly constant wavenumber for all 
values of These horizontal bands of power are the "rings" described by Hill (19SS). 
In polar coordinates, an (Tipse centered at the origin has the form, 

/.• = - (1 - //-) cos-(fJ - 0 q )]~^ (6.1) 

where Aq is half the minor axis (the true wax'enumber); Oq is the angle of the major 
axis (direction toward di^k-centei ): and // is the ratio of minor to major axes (the 
fractional foreshortening). A close inspection of figure 6.3 shows that each of the 
bands exhil^its this oscillator)' belia\'ior of k as a function of 0. 

The \'alue.s of k'o.O^j. and // wei(' fitted to the measured wavenumber k{0) for the 
Pi ring in the freqiKuicy slice // = 3.99 mHz. The function k(0) was defined as the 
wavenumber at peak power in tlie constant frequency slice. It was determined by 
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quadratic iiiteri)olation. 1 he measured fractional foreshortening was, 

fi = 0.956 

This corresponds to a distance from disk-center of 17°, which agrees with an approx- 
imate measurement made' from the full-disk magnetogram in figure 4.2. 

Tlie assumption of a [)lanar atmosphere has been made here, in treating the fore- 
shortening as constant across the field of \dew. If the foreshortening varies appreciably, 
modes of different wa\'enumber will be mixed in the 3-d power spectrum. This will 
smear the trumpet surfaces in (A\//) space in a manner which cannot be corrected. 
For large fields of \'iew, the images must be mapped to "'flattened’' coordinates before 
transforming to Fourier space. The a[)lanarit\' for the data analyzed here was not 
great enough to warrant such treatment. 

The 224" field of view re[)reseiUs 13° along a great circle on the Sun. Since the 
average distance from disk-center was about 17°, plane waves propagating along the 
direction of foreshortening will be mixed o\'er a wavenumber range of ±3.5%. While 
this may seem large, onl\' the waves pro[)agating exactly along a diameter of the disk 
experience mixing to tins extcnit. Wa\'es propagating in the perpendicular direction 
are not mixed at all. The net widening of the ridges is, thus, small. This can be seen 
in figure 6.3. The width of the ridges api)ears independent of 0. The width, therefore, 
is not an artifact of the coordinate' s\'stem. 

6.1.2 Azimuthal Projection 

The power was integrate'd around ellipses in the {kx.ky) plane. To do this, ko was 
computed from equation 6.1 for each ordered pair (kj^^ky) in the 3-d spectrum, and 
power from that bin was added to the' appropriate k bin in the 2-d spectrum. Bins 
in the 2-d spectrum were normalized by (multiplying by) the ratio of the area of an 
elliptical band to the nuniber of bins accumulated from the 3-d spectrum. The area 
A of an elliptical band w<is ap|:)roximated b\*. 



CHAPTER 6. MEASEREMESTS EROM POWER SPECTRA 



1 1 




(Mm.-') 



Figure 6.4: Scaled \'elocit v-power spectrum. A portion ot the collapsed power spec- 
trum lor the lull field ot view, with \’ariable gray-scaling in the k and ly directions. 
Background has been sul>tracted. The noise visible at the top of the spectrum is an 
artifact of the scaling. Frmn Fernandes r/ ai (1992). 

where a and b are the major and minor a.xes ot tlie ellipse; and w is the width of a 
wavenumber bin in the 2-d spectrum. The two-dimensional spectrum P{h,i/) for the 
full-field dopplergrams is ^howii in figure' 6.4. 
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6,2 Parametrizatioii and Least-Squares Fitting 



6.2.1 Fitting Procedure 



The frequencies, widths, and heights of the ridges were measured as functions of 
wavenumber. To do this, constant wavenumber slices of the power spectra were indi- 
vidually fitted by least -scpiares method. Examples of slices at constant wavenumber 
appear in figure 6.5. .A single isolated mode should have a Lorentzian profile. 

Por/2;r 



P(iA) = 



1 ^ 0 ? + rV4 



(6.3) 



where Pq is the mode power: //q is ihc inode frequency; and F is the fwhm. The high 
frequency mock modes dc'scribed by Kumar et al. (1990) should have a sinusoidal 
profile. 



P[iP 






sin”(^\ro/c) 



(6.4) 



r — s\\v{osal c) 

where is the source depth; is the acoustic cutoff; and a is the height of the upper 
reflecting surface. Solar processes such as horizontal flow fields, and measurement 
effects such as seeing and apodization tend to smear the shape to produce a Gaussian 



profile. 

P(u) = (6.5) 

as/'ZTT 

where Pq is the mode power: is the mode frequency; and a = fwhmf2.36. 

The ridges lie on top of a background which decreases monotonically with wave- 
number. In the photosphere, the granulation is caused by overshoot of convective cells 
into the overlying convect ivel\* stable layers. The power, in the photosphere, below 
the /-modes is caused by tins con\*ective o\*ershoot. Extrapolation of calculations of 
Nesis ik: Mattig (19S9) to t ln^ low chromosphere suggests a similar explanation for the 
background power there. 

Because of the uncerl<unty of the physical description for the background, there 
is no proper way to parametrize it. Nevertheless, the background must be accounted 
for to avoid contaminating the ridge fits. It can either be included as an extra term in 
the fits, or fit independently and subtracted from the power spectrum before fitting 
the ridges. 
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Figure 6.5: Constant wavenumber slices of the .3-d power spectrum. Data are plotted 
as points with error bars. Background and individual Gaussian modes are plotted 
as dashed lines; and the sum of background and modes is plotted as a solid line. 
Top left: k = 1.16 Mm“*; Bottom left: k = 5.25 Mm“^; Right (top & bottom): 
k = 2.34 Mm-'. 
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The latter method was chosen l^ecaiise the background could be well characterized 
using the range of frequencies on either side of the ridges. The shape of the background 
Pk suggested a power law parametrization. 



H(,/) = 



( 6 , 6 ) 



where .4^. and Bk are parameters fitted independently for each wavenumber slice. 
The frequency ranges used for the background fits include: (0.12 < ly/mHz < 

1.3) and (T2.0 < ///mHz < 15.0). 

The ridges were fitt(^l after subtracting the background. As can be seen from 
figure 6.5, the ridges o\(‘rlap, and so must be fitted simultaneously. With three 
parameters per Gaussian, and up to eight visil)le ridges, the fits involved up to 24 pa- 
rameters at a time. This large |)arani('ter space coml:)ined with noisy data necessitated 
the de\’elopment of a \*er\' robust fitting algorithm. 

The fitting program is d('scril)ed in the next section. It provides a nonlinear, 
weighted, least-squares fit to ungridded data. The weights used in the least-squares 
fit were the measurement errors. Hence, the fitting procedure minimized The 
error on the power was (estimated from eejuation 2.38. This in turn assumed the 
velocity and intensity errors derived in section 2.1. 

6 . 2.2 Leveiiberg-Marqiiardt Method for Nonlinear Regres- 



The basic algorithm follows the Le\enl)erg-Marquardt procedure (Press et al, 1989). 
The method attempts to minimize the \" in parameter space (a), given N measured 
v^alues ijt ± (Ji at points 



Sion 




(6.7) 



where /y(.r:e) is the fitted function. 



The minimum of the surface is searched for in an iterative manner. Succes- 

sive guesses for the parameter \ ector a are made leased on the gradient and curvature 
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matrix for the surface'. Wlieii the clislance to the minimum, as estimated from 
the curvature, is large, the gradient is weighted more heavily in the guess. When the 
minimum is near, the cui\’ature carries more weight. The algorithm assumes it has 
found the minimum when it can no longer reduce by more than 0.1. The covariance 
matrix for the fit is simply the inverse of the curvature matrix at the minimum. 

To prevent runaways, the parameters are tested after each iteration of the algo- 
rithm in which is reduced. Parameters that are not within set limits are fi.xed at 
appropriate values, and no longer varied in the fit. The remaining parameters are 
reset to their initial valiK's and the procedure starts again from the beginning. 

For Gaussian ridges, l he limits for t he height parameter are t (small, positive) 
and /?,nax (very large). X'aliu's falling outside the limits are fixed at the limit. The 
limits for the width parametc'r are one-half and twice the expected width. Again, 
values falling outside the limits are fixed at the limit. The limits for the frequency 
parameter fall halfway between tlie exj)ected positions for ridges (ao). For the bottom 
(top) ridge in a fit, the lower (upper) limit is defined such that the two limits are 
equidistant from the expected ridge position. \ a lues falling outside the limits are 
fixed at the initial guess uq. 

The reason for this tn'atment of parameter limits is to prevent noisy ridges from 
disturbing the fits of adjacent ridges. Ridges which are too noisy to be properly fitted 
are treated as background. 

For the first spectrum to be fitted, initial guesses at a number of wavenumbers were 
made by e}’e, using cross hairs on a workstation display. Guesses at other wav'enum- 
bers were interpolated or extrapolated from tlie nearest four fitted wavenumber slices. 
For the remaining power spectra, initial guesses were interpolated from the already 
fitted spectra. 

In addition to the l^ackground and ridge fits, the least-squares fitting program was 
used to interpolate and extra [mlate mode frequencies to make the initial guesses, and 
to determine the parameters of eciuation 6.1 in section 6.1.1. 
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Figure 6.6: Power law (\Kponent for background fit. Exponent vs. wavenumber. 

6.3 Results of Fits 

6.3.1 Background 

The exponent to the power law fit is plotted in figure 6.6. The background velocity 
power, integrated over frequency, is plotted in figure 6.7. The spectrum falls off much 
more slowly with wavenumber than the spectrum of convection in the photosphere re- 
ported by Chou et al. (1991 ). The root-mean-square (RMS) velocity is about 820 m/s. 
Nesis Sc Mattig (1989) ha\e measured RMS velocities, due to ‘‘mechanical overshoot”, 
that increase with height in the upper photosphere. Their model predicts the height 
dependence, for line-of-sight velocities, shown in figure 6.8. A smooth extrapolation 
of this curve to the height of formation of the Na Dt line (.500-700 km) yields an 
RMS velocity of 700-1100 m/s, in agreement with our measurements. 




CHAPTER 6. MEASUREMENTS FROM POWER SPECTRA S3 




Figure 6.7; Background velocity power integrated over frequency. Background was 
fitted from the "quiet'’ subfield spectrum. From Fernandes et al. (1992). 

6.3.2 Mode Frequencies 

Frequencies of ridges from the full-field velocity and intensity spectra, and the quiet 
subfield velocity spectrum are tabulated with their respective statistical errors in 
appendi.x A. The velocity points are plotted in figure 6.9. All tabulated points 
correspond to fits in which the offset, width, and height of the Gaussian converged 
to physically reasonable values, and the integrated power under the Gaussian was at 
least 4 standard deviations above the background. 

All six fitted p-mode ridges extend up to at least 7.0 mHz. The ps and pe ridges 
extend up to a frequency of 9.-5 mHz. There is no deviation from a smooth extrap- 
olation of the ridges at lower frequencies. Libbrecht, Woodard & Kaufman (1990) 
have made measurements that overlap with ours at frequencies up to 5.3 mHz. As 
can be seen from figure 6.10. the fitted frequencies are consistently higher than those 
of Libbrecht et al. by an average of (35 ± 13) pHz. This difference, however, is only 
a quarter of the frequency bin width, and is indicative of the systematic error in the 
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Figure 6.8: Calculated Doppler velocity vs. height for the Fe I (5123.75 A) line. From 
Nesis <Sc Mattig (1989). 

measurement. 

The /-ridge is visible at wavenumbers up to 5.57 .VIm“b There is, however, a 
systematic discrepancy between frecpiencies fitted from the velocity and intensity 
spectra as can be seen from figure 6.11. Frequencies measured from the velocity 
spectra are consistently higher than those from the iiitensity spectra by up to 50 /iHz 
at wavenumbers below 2.5 Mm~’. This difference increases to ~ 100 /zHz at 4 Mm~^ 
The discrepancy may in part be due to a difference between background shapes in 
the velocity and intensity spectra. 

Libbrecht et al. (1990) report a Gaussian background, centered near 3 mHz, in 
their velocity spectrum. Inclusion of such a background in our velocity fit would bring 
the /-ridge frequencies into better agreement with the intensity data at wavenumbers 
below 1.1 Mm“\ but would push them farther apart at higher wavenumbers. Fitting 
such a background function, moreover, would require a longer time series with corre- 
spondingly better signal-to-iioise ratio. It is noteworthy that the /-mode frequencies 
measured from the intensity spectra differ with those reported by Libbrecht et al. by 
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Figure 6.9: Measured mode frequencies (from velocity spectra): i/ vs. logj^ Fre- 
quencies of the / and p\-]M ridges are plotted with 1 a error bars. For k < 3.82 Mm“\ 
frequencies are taken from the full-field spectrum. At higher k, values from the quiet 
subfield are plotted. 
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Figure 6.10: Comparison of measured /j-mode frequencies with Libbrecht et al. (1990). 
Frequency difference vs. k is plotted for pi,P 2 ,P 3 (from top down). 
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Figure 6.11: Comparison of /-mode frequencies measured from velocity and intensity. 
Frequency difference vs. h. 

under 30 //Hz. as shown in figure 6.12. 

The weighting of magnetic regions is different for intensity and velocity spectra. 
However, as shown in section 6.3.5, the discrepancy is not caused by a difference in 
mode frequencies in the magnetic and nonmagnetic regions within the field of view. 
Regardless of which set of /-mode frequencies we use, the measured frequencies lie 
below the theoretical frequencies given by equation 1.9 at wavenumbers greater than 
1 Mm~'. The difference is 100-200 //Hz at 4 Mm“’. The fractional difference between 
measured and theoretical frequencies is plotted in figure 6.13. 

6.3.3 Mode Widths 

Ridge widths are plotted as error bars to the frequencies in figure 6.14. The ridge 
widths at low frequency are due to the apodizing function. At high frequency, the 
ridges are broader than expected for mock modes. The meaisured widths probably do 
not represent the true widths. The width is the most sensitive parameter to noise in 
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Figure 6.12: Comparison of /-mode frequencies from full-field intensity spectrum with 
Libbrecht et al. (1990). 

the fitted power spectrum. 



6.3.4 Mode Power 



Mode power is defined as the integral under the fitted Gaussian. 



2tt 



( 6 . 8 ) 



where h is the height of the ridge; and w, the fwhm. The error on the power is given 

by, 

(6.9) 



(Tp = P] 



+ 



+ 2 






h J \ W J hiv 
where <Thw^ is the covariance of the height and width. The mode power is plotted 
as a function of frequency in figure 6.15. The p-mode power has the same frequency 
dependence for all p-ridges. There is very little explicit wavenumber dependence. 
Power increases as i/''^ to a maximum at .3.4 mHz, at which point it turns over and 
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Figure 6.13: Comparison of measured /-mode frequencies with theory. Frequencies 
from the full field intensity spectrum are compared with the relation = gk {g = 
274.1 m/s^). 
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: Frequency and ridge width vs. k. Top plot shows /-modes and even 
. Bottom plot shows odd ordered p-modes. Error bars represent the 
! fit. 
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falls off again as The /-modes have less power at a given frequency than the 

p-modes. As surface gravity waves with very little compression they are physically 
different from the p-modes and are not expected to show the same power spectrum. 

An MTF correction was made to the power plotted in figures 6.7 and 6.15 by com- 
bining the calculated MTF of the telescope and CCD with the estimated atmospheric 
MTF. See sections 3.4 and 5.4.1. 

6.3.5 Comparison of Magnetic and Nonmagnetic Regions 

Mode frequencies were compared for the quiet, canopy, and sunspot subareas shown 
in figure 6.1. The frequency differences are plotted in figures 6.16 and 6.17. The 
frequencies of ridges in . each subarea agree with those of the other subareas to within 
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Figure 6.16: Mode frequency differences for canopy and quiet subfields. Top left, 
/-mode; top right, pi; bottom left, pY, bottom right, p 3 . 

the scatter of the data (±40 /tHz). 

The experimental scatter is much larger than the error bars derived from the fit. 
This may simply indicate an underestimation of the velocity error, or it may be a 
result of solar noise. Modes of finite lifetime have an intrinsic frequency uncertainty. 
Determination of the frequency of the cavity resonance j/q requires observations for 
several mode lifetimes. The 2.4 hr observing run probably constitutes less than one 
mode lifetime over much of the wavenumber range. Hence, the measured frequencies 
will be distributed about uq in a manner which depends on the mode lifetime, and 
shape of the wavepackets. 
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Figure 6.17: Mode frequency differences for sunspot and quiet subfields. Top left, 
/-mode; top right, pi; bottom left, p 2 \ bottom right, ps. 



Chapter 7 

Spatial Distribution of Power 



The spatial clistril)ution of power was measured for groups of modes, as well as for the 
sub-fundamental region of tlu^ power s])eetrum. To do this, the set of dopplergrams 
V{x,t) was masked using 11 (T) from equation 2.33, and a 3-d FFT, calculated for the 
full field of \'iew. A spectral mask was then applied to the Fourier transform 

to remove all but a small range of wavenumbers of a given ridge. An inverse Fourier 
transform of this masked 3-d spectrum yielded velocity for a specific set of 

modes. 

\^,(7./) = ^-q.U(C;.)^qir(f)b(.r,0}} (7.1) 

The mean power was cal< ulated at each pi.xel to make a \’elocity power map Pm{x). 



I Jo 
I 'J-l'uvn 



du 



(7.2) 

(7.3) 



where denotes the spatial, inverse Fourier transform. 



7.1 Mode Masks 

Masks were constructed fur the / through /> 3 - modes in 6 wavenumber ranges. In addi- 
tion, two sub-fundamental and a liigli frequency set of masks were made for the same 
wavenumber ranges. The mask limits are shown in figure 7.1. Three-dimensional 
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masks were made from these 2-d limits by rotating the loci in the azimuthal direc- 
tion. To account for the foreshortening, the rotation followed ellipses described by 
equation G.l. ! 

I 

7.2 Results 

Two sample power maps are compared with a magnetogram and an H a image in i 
figure 7.2. Only the central disks oi the frames are visible due to the apodizing 
window. The subi-/and /-mode maps W(u*e both made for the wavenumber range 
1.95 < A’/Mm”* < 2.51. I his is (Hpii\‘alent to a wavelength of about 4", or 3 Mm. 
The mottled appearance of the power maps is an artifact of random interference of ’ 
waves in a finite time scnies. The subi-/map shows little variation other than an | 
overall trend of higher i)ower in the upper right hand corner. The /-mode map shows 
a depression in power in the active region. 

A selection of power maps is shown in figure 7.3. An estimate of the RMS velocity ; 
in each map can be obtaiin^d from figure' 7.4. The power in high wavenumber modes 
is depressed in the active* re'gion. At low wavenumbers, the active region has little 
effect on the i)Ower. 14 k' subi-/and siibo'/maps show no correlation of power with ' 
magnetic field at an\’ wa\ eiumibc'r. 44ie high frequency maps, though not shown in i 
figure 7.3, looked very similar to p-mode maps of the same wavenumber range at 
lower frequenc}'. There is little difference between the /, pi, p 2 , and pa maps for a 1 
given wavenumber range. 

A horizontal slice through the /-mode map from figure 7.2 was made above the 
sunspot. The slice was low |)ass filtered using a Gaussian impulse response function 
(IG bin fwhrn). See figiiie 7.5. I'his low pass filtered slice is compared with the 

I 

magnetogram signal in figure* 7.G. Idie power depression is significantly broader than 
the spatial extent of the plioticspheric magnetic field. It is narrower than the canopy 
field at ~ 1200 km, as can lx* s(*(*n from the II a image in figure 7.2. The filaments 
visible in 11 n are streams of gas traveling along magnetic field lines. The filaments, 
thus, trace out the field lines which form a canopy above the photosphere and lower 
chromosphere. 
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k (Mm 

Figure 7.1: Masks for power distribution maps. Measured mode frequencies are 
shown as dots. Mask limits are shown as solid lines. The high frequency (HF) masks 
covered the range 6.0 < u/mHz < 10.0; the p-mode mask limits were set halfway 
betw'een modes; the /-mode limits were set equidistant from the /-mode; the sub 2 -/ 
band e.xtends from the /-mode down to = 0; the subi*/band is the same width as 
the /-mode band. It is a subset of the sub 2 -/band. 
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Figure 7.2: ( oinparisoii ul maps with magnetic Helds. Top left: log power 

lor si\h\-J (1.9-) < /'/Mm ’ < 2..51). lup right: absolute \’alue of magnetogram (log 
scale). l3oltom h*lt: 11 n — lOU mA image. Bottom right: log power tor /-mode 
(1.95 < /'/Mm ^ < 2..'M). 1 1 k‘ small tick marks are 
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Figure 7.3: Power (list ribiil ioii lor selected ina.sks. Log |)ower is dis|)layed for sub-j-/, 
/> Pi- Pi- I>3 niajis with wiuenmiiLer range.s: 

[1.95.2.31], [I.IT. 1.5(ij. [0.78.1.17]. [0.51 . 0.78] .Mm"'. 
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Figure 7.4: RMS velocity vs. k. Curves fitted from quiet subfield velocity spectrum. 
Curve labeled sub-/ is taken from the background fit. 

It is surprising that the sub-fundamental power is not correlated with magnetic 
field. The RMS convective power in the photosphere is a factor of 2 lower in plage 
(network field) than in nearby quiet regions (Title et al. 1992). Moreover, the sub 2 -/ 
band lies just below the /-mode band, but shows completely different behavior. At 
the same wavenumber, and almost the same frequency, oscillation amplitudes are 
reduced in active regions, whereas non-oscillatory motion is not. 

It is impossible to distinguish between instrumental backgrounds and solar back- 
grounds from a single data set. It is difficult, however, to come up with an obser- 
vational effect that would produce the measured background. Filter transmission 
variation should affect only the very low wavenumber part of the spectrum. Seeing 
variations should be independent of frequency; whereas, the observed background 
decreases with frequency. 
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Figure 7.5; Slice of /-mode map. Horizontal slice from the [1.95,2.34] Mm"‘ map. 
Also shown is a low pass filter of the slice. 
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Figure 7.6; Comparison of /-mode slice with magnetic field. Dashed line shows low 
pass filter of slice (arbitrary scale). Solid line shows magnetogram signal (arbitrary 
scale). 



Chapter 8 

Conclusions and Discussion 



8.1 Surfaces and Cavities 

The existence of a chroniuspheric cavity is expected to produce a resonance consist- 
ing of a liorizontal band of pcjwcr and "avoided crossings'’ at about 5.5 mHz. (See 
figure 1.4.) Tfiis ca\*ity is recpiired l:>y the model of Balmforth Sz Gough (1990) to 
explain tlie existence of high freciuenc}' inodes. The /;-mode and /-mode ridges shown 
in figures 6.4 and 6.9 are ( learly visible at 5.5 mHz, but there is not a horizontal band 
of power nor are there "ax oided crossings". The clear absence of this resonance in our 
data shows that there is no chromospheric caxhty. Either the reflection coefficient at 
the base of the corona is loo low or the surface is too irregular to form a good cavity. 

This result supports the alternate hypothesis of Kumar et al. (1990), that the 
high frequency ridges arc' a piirelx' geometrical artifact. The result is in agreement 
with the conclusion of Woodard Libbrecht (1991); and it contradicts the claim of 
Dame et al. (1984) that there is a horizontal ridge of power at 5.5 mHz. Our power 
spectrum has a significantly better signal to noise ratio than that of Dame et al. 

The /-modes are sensitix'e lo conditions in tlie chromosphere and transition re- 
gion. The fact that /-mode freciuencies are lowered at high wavenumber may indicate 
interaction with cliromospheric magnetic fields, energy leakage into the corona, or 
horizontal inliomogeneitic's in tlie chromosphere-corona boundary. The last of these 
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three possibilities is iiUcK'sting from the standpoint of acoustic cavities. Perhaps ir- 
regularities in the interfac e l^etween the chromosphere and corona are responsible for 
both lowering the /-mode fK^iuencies, and for reducing coherent reflection of acoustic 
waves. 



8.2 Waves in Magnetic Fields 

There is no frecpiency shift for / and p-modes in magnetic network fields of 100 — 
— 1000 gauss to within lh(‘ scatter of tlie data 50 //Hz). This limit is not stringent 
enough to determine' wlu'the'r the frequency shifts measured by Woodard et al. (1991) ; 

I 

result from small shifts over the whole disk, or larger local shifts, in active regions, i 

I 

which contaminate the full-disk measurement. The modes measured by Woodard et 
al. have widths of order 1 //llz. Hence, to contaminate such a measurement, mode 
frecjuencies in actix'e regiuus would have to be distributed within a few /^Hz of the 
quiet sun values. A measiii'eiiK'iit of such precision would require observations of at • 
least 20 hr duration, assuming the wa\*es are coherent for the observing period. At 
La Palma, 20 hours of obser\-atioiis can be made in two days. Hence, the waves would 
add coherently if the mode lifetiiiK' is > 31 hr. If the mode lifetime is less than this, 
the signal to nois(' ratio iucK'ast^s only as the square root of time, so a much longer 
observing time would be n(‘C('ssar\'. 

Power maps for oscillations in the low chromosphere show that, at high wavenum- j 

I 

ber, wave amplitudes ar(' su|)pressed 15 Mm outside of the photospheric magnetic 
field. (See figure 7.6.) In (he photosphere, the power is suppressed within one (hori- 
zontal) wavelength from the held (Tarbell ci al. 19SS). The vertical distance between 
the photospheric and chromospheric observations is under 0.5 Mm. and the (true) 
wavelength of the oscillations is -- 2 Mm. Ditfraction cannot account for the spread- 
ing of the region of the powt'r depression because the waves have traveled less than 
one wavelength in the inl('r\'al o\-er which the spreading occurs. The power suppres- 
sion must, then'fore, be produccxl by interaction of the oscillations with the field in 
the layer in where the wart.'> are obserred. 

The area containing magnetic flux was not measured directly at the height of the 
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observations. Its extent ran be inferrecb however, from the niagnetogram made in 
the mid-photosphere, ami tlie II a image in the mid-chromosphere. The field lines 
spread out rapidly abov(' the photosphere. Hence, the extent of the field at the 
observed height should h(' intcu'mediate between the limits of the mid-photospheric 
and mid-chromospheric fitdds. This is consistent with the above conclusion. 

The /-modes propagate without compression, so they produce no pressure fluctu- 
ations in the region where oscillations are suppressed. The mechanism suppressing 
the wave motion must, therefore, act through the velocity' fluctuations. 



8.3 Directions for Future Work 

The frequency ineasurenK'iits presented here show the existence of high frequency and 
high wavenumber oscillations, bhey show general characteristics of the dispersion 
relations such as the smooth continuation abo\*e the acoustic cutoff. They are not, 
however, high precision measurcmients in comparison with existing low wavenumber 
observations. Precision of the measurements can be improved with longer observing 
times. This would impro\e the signal to noise ratio, and enable better fitting of the 
background and ridges. With better frec|uency resolution, it would become necessary 
to correct the fitted fre(|uencies for their weighting with mode power (Korzennik 
1991). 

Important characterist ics of the solar oscillations can be inferred from their phases 
(Deubner Sz Fleck I9S9, 1990; Fleck Sc Deubner 19S9). Trapped and nontrapped 
oscillations can be distingiiislied by the phase clilTerence between velocity and intensity 
measurements. Propagating and evanescent waves can be distinguished by the phase 
lag between velocities at dirferent heights in the atmosphere. An attempt was made 
to measure velocity-intensit}’ pha.se differences. The attempt failed, possibly because 
of the large height range included in the measurement (see section 1.4), and possibly 
because of the large number of modes included in a single frequency bin (Braun et 
ai 1992a). 

The interaction of magnetic fields with oscillations can be better quantified by 
measuring both wa\'e amplitudes and magnetic fields at a range of heights in the 
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atmosphere. Improvemeiils that ha\'e been made to the tunable filter system since 
the 1990 ol^serving season will allow Do[>pler measurements restricted to either the 
core or wings of the Na Di line. In addition, magnetograms can be made in the core 
or wings of the line b\- finding the difference between dopplergrams made with left 
and right circularly polarizc'd light. 

Braun ct at. (1992b) report halos of enhanced high frequency acoustic power 
around active regions. A search was made for such halos in power maps made with 
the HF masks shown in figure 7.1. Xo enhancement of power was found. This does 
not contradict the result of Braun ct at., as they measured intensity fluctuations in 
the photosphere, in a lower waveiuimljer regime. It would be interesting to repeat 
the analysis with continuum brightness measurements for a field of view of at least 
6 arc minutes. 

The cause of the sub-fundamental background is, currently, unknown. Experi- 
ments should be doin' to discriminate Ix'tween observational and solar backgrounds. 
For example, a set of filtergiams taken at fixed wavelength with a high cadence 
could be used to make artificial dopplergrams, and check the effect of seeing on back- 
ground. .Another set of filiergrains taken at the nominal, observing wavelengths could 
be made, using integrated sunlight, to check the effect of filter- transmission variation 
on background. 

An attempt was made to find the oscillation sources described by Brown (1991), A 
temporal!}’ high-pass filtc'red movie of velocities was scanned visually for evidence of 
spatially and temporall}’ localized sources of power. No good candidates were found. 
A better method might b(' to construct the two-dimensional acoustic Pointing vector 
described by Braun ct oL (1992b). Large positive divergences of this vector field for 
high-pass filtered data could [)oint out oscillation sources. 



Appendix A 
Frequency Tables 



Table A.l: Mode frequency u ±(t„ (/‘Hz) as a function of wavenumber k (Mm”^) — 
measured from dopplergrams. For k < 3.S2 Mm~\ frequencies are taken from the 
full-field spectrum. .-\t higher k. values from the quiet-subfield are tabulated. 
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Table A. 2: /-mode frequency j/ ± cr^ (/^Hz) as a function of wavenumber k (Mm *) 
— measured from intensities. 
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